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The possibility of the extension of spatial diffeomorphisms to a larger family of symmetries 
in a class of classical field theories is studied. The generator of the additional local symmetry 
contains a quadratic kinetic term and a potential term which can be a general (not necessarily 
local) functional of the metric. From the perspective of the foundation of Einstein's gravity 
our results are positive: The extended constraint algebra is either that of Einstein's gravity, or 
', ultralocal gravity. If our goal is a simple modification of Einstein's gravity that for example 

makes it perturbatively renormalizable, as has recently been suggested, then our results show 
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that there is no such theory within this class. 
I. Introduction 



What makes an n + 1 dimensional (pseudo) Riemannian manifold the relevant structure in the 
description of the time evolution of Riemannian ^-geometries? In the Hamiltonian formalism, 
which is based on a spacelike foliation of the spacetime, the invariance of Einstein's gravity under 
nonspatial diffeomorphisms - diffeomorphisms that cannot be restricted to the spatial slices of 
the foliation - is somewhat hidden, whereas the spatial diffeomorphism symmetry is as manifest 
Q\ [ as in the Lagrangian formalism. More concretely, it requires some calculation to realize that the 
infinitesimal symmetries that the Hamiltonian constraint generates on the solutions to the equations 
■^j- ■ of motion correspond to infinitesimal deformations of the foliation. Since the relationship between 
(sj . the transformations generated by the constraints and the diffeomorphisms of the underlying n + 
1 dimensional manifold is not straightforward, the question naturally arises if this relationship 
is necessary at all. To put it differently: Do spatial diffeomorphisms allow for a further local 
symmetry so that the symmetry algebra is different from that of general relativity or ultralocal 
gravity? We insist on general covariance, however, we use this term in a weaker sense than usual: 
The spatial metric and the conjugate momentum are the only canonical variables, and the local 
constraints are constructed out of them in a form invariant way. 

This is not the first time that the possibility of reducing the set of postulates that lead to Ein- 
stein's gravity is considered. 1 ^ However, we are interested in this question not only from the 
perspective of the foundation of Einstein's gravity. Another motivation is a recent proposal by 
Hofava for a modified theory of gravity in which Lorentz invariance is given up in the quest for 
an improved UV behavior. 5 Although this model does not seem to be physically relevant, it puts 
our question in a new light: Perhaps it is possible to reduce the postulates of general relativity so 
that they are realized not only by general relativity and ultralocal gravity, but also by some other 
theories which might have the properties that motivated Hofava's modification. First we briefly 
describe Hofava's proposal, and set up the framework in which we can look for theories that can 
be relevant for gravity. 
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The starting point of Hofava's construction is the action of general relativity in terms of the 
ADM Lagrangian: 

Sadm= / ^adm= / Nyfg(R + KyK ij -K 2 ), Ky = ±-(gy - 2V (l -N /) ), (1) 

Rxl Rxl 2JN 

where I is the space, and R is the time. The field configuration on the spacetime manifold is given 
by (g(/,N^), W = (N, N ! ), where N is the lapse, N 1 are the shift functions, and gy is the spatial 
metric. In (1) R is the Ricci scalar of gy, and V is the covariant derivative on L compatible with 
the spatial metric. Ky and K are the extrinsic curvature and its trace, and gy is the time derivative 
of gy. Latin indices denote spatial components. A configuration in the phase space V is given by 
(gy, Ky). To avoid indices, we will denote the elements of V by (g, n), where n is the momentum 
conjugate to g. 

Spatial diffeomorphisms are still local symmetries in Hofava's proposal as in Einstein's gravity. 
What remains from the nonspatial diffeomorphisms is the invariance of the action under space in- 
dependent, i.e. foliation preserving transformations. For further restrictions on the possible action, 
Hofava first appeals to an effective field theory argument: Under an anisotropic scaling symmetry, 
the theory has a UV fixed point, and apart from the UV potential the Lagrangian contains only 
relevant couplings under the scaling symmetry. For pragmatic reasons, a further principle is intro- 
duced. The theory should satisfy a "detailed balance" condition, which restricts the form of the 
potential term in the action: 

S= f*ANV?(^« H -^„g. (2) 

where W and the tensor G iJkl are some functionals of the spatial metric, and k is a constant. The UV 
potential is proportional to C ij Cy, where C ij is the Cotton tensor. The detailed balance condition 
limits the number of the possible relevant couplings under the anisotropic scaling symmetry if we 
assume that this condition is preserved by the renormalization group flow. The action suggested 
by Hofava is (2) with 

G W =g * g jl + g U g jk_ Ag u g H t (3) 

W = ^ J d 3 x Tr (r A dr + A r A rj + // J d 3 x yfg(R - 2A W ), (4) 

where A, w, /u, and A w are constants, and the first term in W is the gravitational Chern-Simons 
term, expressed in terms of the Christoffel symbols 7"}* of gy. The functional derivative of this term 
with respect to the metric is proportional to C ij . 

By discarding possible spatial boundary terms, (2) and its later generalizations can be got from 
the following action given in the Hamiltonian formalism: 

S = fdtf(g-7r-N»J^[g,n]), (5) 

where at any time = (J^, are functionals of the spatial configurations of (g, n). This action 
is defined on a class of time dependent configurations (t, x) i-» (g(t, x), n(t, x), N M (t, x)), t e (t\, ti), 
x e L. Any field that can be an instantaneous configuration in a function within this class will be 
referred to as kinematically possible. A field that can occur as an instantaneous configuration in a 
(local) solution to the equations of motion of (5) will be called dynamically possible. Let V be the 
constraint submanifold of F, which, by definition, consists of the dynamically possible canonical 
variables. We will assume that the infinitesimal transformations generated by first class constraints 
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on f can be integrated to a symmetry group CAN. Subscript ADM will denote the corresponding 
objects in general relativity. 

In the nonprojectable version of Hofava's theory N can be any smooth function apart from 
some conditions on its asymptotic behavior. In this version the equations of motion include the 
constraints = on V, which are not first class. 6 We mention three different approaches to this 
problem: 

(i) We can get first class constraints by imposing further constraints on P This can be implemented 
by adding constraint functions with multipliers to (5). The new symmetry group CAN' is an 
extension of the spatial diffeomorphisms by local symmetries. The new constraint submanifold 
P £ r is a submanifold of V, so P and especially the manifold P/CAN' of orbits are too small 
compared to these objects in general relativity. Thus this option cannot be relevant for gravity. 

(ii) The nonspatial diffeomorphisms in Einstein's theory are symmetries of (2) only if they pre- 
serve the foliation. These symmetries can be implemented on the action if the set of kinematically 
possible configurations is invariant under such transformations. For example N can be kinemati- 
cally possible only if it is a constant function, and thus the spacetime configurations of N depend 
only on time. This is also called the projectable version of the theory. Here the constraints are first 
class since apart from the constraints implied by the spatial diffeomorphism symmetry, there is one 
single global constraint which is invariant under the transformations of the canonical variables by 
spatial diffeomorphisms. 

(iii) It is not necessarily a pathology if the constraints are not first class. It just means that (g, n, N^) 
is not dynamically possible for any kinematically possible N and (g, n) e P It is possible that we 
do not have to modify the theory as in (i) if our only goal is that the local constraints = 
are preserved by a nontrivial time evolution, for which all the constraints need to be satisfied by 
a nonzero N for any (g, n) e P A dynamically possible configuration of N[g Q , 7r ] for a given 
(g , 7r ) G F satisfies { J^[g, n](x) , f L N[g , n ]Jif[g, n] }\ g=go ^ =jro = 0, where the Poisson bracket 
is taken with respect to (g, 7r). On the other hand, it is necessary that N[g , 7r ] satisfies this equa- 
tion in order for SF[g Q ,n ] = {F[g,n] , f L N[g Q ,7ro]Jf[g,7i], }\ g=gQ ^ =7TQ to define an infinitesimal 
symmetry on P It was argued in [7] that in a class of theories which includes Hofava's proposal 
the only solution to this equation is N[g , tt ] = for generic coupling constants and (g , n ) e P 
Thus there is no symmetry associated with Jff. 

There are several reasons why the lack of the local constraint J$? = or the corresponding 
local symmetry is undesirable. If A £ | in (3) then the local scale factor of the spatial metric is 
dynamical, i.e., using the variables <p := | lng and gy := gye'^ instead of gy in (2), we can see that 
KjjG ijkl K k i is independent of <p if A = |, whereas it is quadratic in </> for any other A, and the acceler- 
ation of (j) appears in the field equations. If A > | in (3), then the signature of G ijkl as a metric on the 
symmetric rank two tensors is (-, +, +, +, +, +), and the term in (2) proportional to 4> 2 is negative. 
Such a mode usually leads to the loss of unitarity in quantum theory unless it is unphysical, i.e., 
there are gauge symmetries, and there is such a gauge condition that fixes this mode apart from 
some nonlocal degrees of freedom, which is not the case if the only local symmetries are the spa- 
tial diffeomorphisms. In Einstein's gravity it is the space dependent nonspatial diffeomorphisms, 
which are generated by the Hamiltonian constraint, that make the gauge choice K\ = admissible 
at least for spatial metrics that satisfy some asymptotic conditions in suitably chosen coordinates 
(Dirac's maximal slicing gauge, see for example [8]). Note that A changes between A = | and 
A = 1 if general relativity is to be recovered from (2) in some limit, and the local scale factor of 
the metric is a physical mode with negative kinetic energy for such A if there is no symmetry that 
replaces the temporal diffeomorphisms. Even if it was possible to define a unitary quantum theory 
with such a physical mode, its relationship with Einstein's gravity, where this mode is unphysical, 
would be unclear. 
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The lack of a local symmetry that replaces space dependent temporal diffeomorphisms also 
makes the physical interpretation of certain solutions problematic. Let /i = in the potential (4), 
and let us break the detailed balance condition, which seems to be necessary in order that general 
relativity be some limit of the modified theory. 9 If we add the Ricci scalar of the spatial metric to 
(2) so that the limit w —> oo yields the ADM action (1), we get 

S = j dtd'xNyfg^KijG^'K^ - ^C ij Qj + jR^, (6) 

where we used that C ij is traceless, so Gy k iC ij C kl = 2C ij Cjj. A spherically symmetric spatial met- 
ric is conformally flat, and the Cotton tensor is zero for such metrics in three dimensions. For 
(gij, N, 0), where gy is conformally flat and gy = 0, the field equations of (6) reduce to the same 
equations as those of the ADM action (1). So the configuration (gy,N,N) obtained from the 
Schwarzschild spacetime metric in the Schwarzschild coordinates is a solution to the field equa- 
tions of (6). The horizon is a coordinate singularity in Einstein's gravity, but here it appears to 
be a physical singularity, where the spatial metric changes its signature. In general relativity the 
coordinate singularities can be removed locally by taking a family of timelike geodesies heading 
towards the singularity, and using their proper time as the new time coordinate. Such a family 
of timelike geodesies can be arbitrarily approximated by physically realizable observers (by some 
matter distribution). In Hofava's theory, with spatial diffeomorphisms as the only local symmetries, 
there is a preferred foliation of the spacetime. The label of the spacelike hypersurfaces serves as a 
time coordinate. If this time coordinate has any physical meaning in the sense that there are phys- 
ically realizable observers whose time arbitrarily can approximate it, the spatial metric changes its 
signature at some points according to such observers. This phenomenon is hard to interpret. 

Of course it is possible that such observers are not realizable. For example the acceleration 
of a stationary particle may be unbounded as its radial coordinate approaches the location of the 
horizon. Note that even if this is true in Einstein's gravity, it does not necessarily hold in Hofava's 
modification. It is not obvious how Hofava's proposal incorporates matter, and if test particles 
respect the geodesic principle. The geodesic principle is true if the tensor whose nonzero value 
indicates the presence of matter is obtained by varying a diffeomorphism invariant action with 
respect to the metric, and it satisfies the dominant energy condition. (See [10] for the details, in- 
cluding the precise meaning of the geodesic principle here.) Even if we reconstruct a spacetime 
metric from the field variables in Hofava's theory, this result does not apply. Nevertheless, assume 
that the geodesic principle holds. It is a natural requirement that there be a time coordinate which 
can be measured by a procedure which is physically admissible in the entire space that should con- 
tain the region where the spatial metric is properly Riemannian. So we are in the puzzling situation 
that the acceleration of stationary test particles can be arbitrarily high in any time coordinate that 
is defined by physically feasible instructions. 

As mentioned earlier, Hofava proposed a UV theory that has a local scale invariance. The 
action is (2) with only the least relevant potential term: 

5 UV = j d ? d\NVi||(^'-^ 2 )-^C%j, (7) 

which is invariant under (gy, N, N) — » (gye 2co ,Ne 3< °,N l ), where co is an arbitrary smooth function 
on L. The Cotton tensor transforms as C ij —* C i -'e~ 5(0 . In order to implement this symmetry, 
the space of the kinematically possible configurations of N should be invariant under the scale 
transformation. This space is big enough to get a local constraint by varying the action with 
respect to N. The constraints on the phase space are not first class (see [6] for direct computations 
for the special case when the action is the same as (7) except for the kinetic term, which is taken to 
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be the same as in the ADM action (1)). It is a relevant question what generalizations of the initial 
class of theories might allow for a first class constraints, or if it is possible to realize the local scale 
transformation so that this additional local constraint is absent. 

We have argued that the plausible candidates for a modification of Einstein's gravity seem to 
be the ones that admit an extension of the spatial diffeomorphisms with a further local symmetry. 
Actually, this could be the guiding principle in our construction of theories in which only the 
spatial diffeomorphisms are left from the symmetries of general relativity: We require that the 
temporal diffeomorphisms are replaced with a new local symmetry. We have shown above some 
physical motivations for this requirement, and we think it should replace the detailed balance 
condition whose status in Hofava's original proposal is dubious. It was introduced only to avoid 
the proliferation of the possible terms in the Lagrangian, and in fact it turns out to be too restrictive 
if the theories that respect it are designed to reproduce Einstein's gravity in some limit. Breaking 
the detailed balance condition in order to overcome the latter difficulty raises the question that the 
imposition of this condition initially intended to answer: Why do we not have to switch on the 
other relevant terms that break the detailed balance condition? There are intimidatingly many such 
terms. 

Once we accept the need for extensions of spatial diffeomorphisms, we have to choose a frame- 
work in which they are easy to find. The symmetries of Einstein's gravity take a simple form in the 
Lagrangian formalism. They are diffeomorphisms of a four-dimensional manifold. In the canon- 
ical formalism this simplicity is lost. The generators of the infinitesimal transformations do not 
even form a Lie-algebra, owing to the field dependent structure functions of their Poisson algebra. 
Nevertheless, nothing is lost from the local symmetries of solutions in the Lagrangian formula- 
tion as we pass to the Hamiltonian formalism in the following sense. Consider & — > r ADM , 
where the bundle manifold & consists of maximal globally hyperbolic vacuum spacetimes. r ADM 
and CAN adm are as before. See [11] for more details about the construction of these manifolds. 
Let Diff be the group of not necessarily metric independent spacetime diffeomorphisms. Then 
II : #/Diff -» r ADM /CAN ADM , TT(Diff • g) := CAN adm ■ <9) is a well-defined bijection. 11 ' 12 If 
we accept that plausible theories share this property with general relativity, the canonical formal- 
ism is an appropriate framework for our investigations, and we can use it without loss of generality. 
In our case the possible extensions of the spatial diffeomorphism algebra are yet unknown, so we 
do not know a priori in what formalism, if any, the symmetry group takes such a simple form as 
the spacetime diffeomorphisms in the Lagrangian formalism of general relativity. 

Similarly to Hojman, Kuchaf, and Teitelboim, 1 3 we assume that the canonical variables are the 
spatial metric and its conjugate momentum, the same as in Einstein's gravity. The basic difference 
between their analysis and ours is that they assumed that the symmetry algebra is the canonical 
representation of the surface deformation algebra, which is the algebra of deformations of an n di- 
mensional spacelike surface embedded in an n + 1 dimensional manifold, where the deformations 
are induced by the infinitesimal diffeomorphisms of the n + 1 dimensional space. In their case the 
Poisson algebra of the generators was completely known. In our case the n dimensional diffeomor- 
phisms of the surface are part of the symmetry transformations, but unlike Hojman, Kuchaf, and 
Teitelboim, we make no assumption on the form of the entire symmetry algebra. We are interested 
in the existence of algebras different from the representation of the surface deformation algebra 
rather than the uniqueness of a realization of a fully specified algebra, which is the subject of the 
analysis in [1-3]. 

Section II and III contain some preparatory definitions and a discussion of the momentum 
constraint. Section IV analyzes the possibility of extending the supermomentum, which generates 
the spatial diffeomorphisms, by an additional generator function, whose form is taken to be a 
local kinetic term quadratic in the momentum n ij together with a potential of quite general form. 
Under some simplifying assumptions on the structure functions, we find no extension of the spatial 
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diffeomorphisms within this class, apart from the symmetry algebra of general relativity and its 
ultralocal truncation. Hofava's proposal belongs to this class. 

In Section V we modify Hofava's candidate for the UV fixed point theory. The phase space of 
this modification is reduced to the variables which are invariant under the local scale transforma- 
tion. The constraints on the the phase space are first class, and they all generate symmetries on the 
solutions. In Section VI we mention some possibilities that have not been ruled out by the negative 
results of Section IV. 

Several appendices are attached to the paper. One of their purposes is to maintain the mathe- 
matical rigor of our analysis. Appendix A is devoted to some technical assumptions, which could 
be replaced by stronger locality conditions on the quantities arising in the Poisson algebra of the 
constraint functions. Nevertheless, this appendix helps us see what properties are really impor- 
tant for the final conclusion, and we can keep closer to the generality of [1], where no locality 
assumption was made on the metric dependence of the elements of the Poisson algebra. In light of 
speculations that events might not play any essential role in the theory of gravity, the possibility of 
nonlocal gravitational potential that partially smears out events is worth considering. Appendices 
B, C, and D describe some mathematical properties of the constraint functions, which the reader 
might find interesting in their own right, and which could be used in a further investigation of the 
symmetry algebra. 

II. Technical preliminaries and notations 

Throughout the paper a symbol with an arrow on it denotes a vector field, a boldface letter stands 
for a general tensor density, or sometimes in the appendices, a tensor field valued linear map. 
Normalface letter with indices denotes the components of a tensor density or the tensor density 
itself, without indices, its trace (e.g. n, n ij , n). 

In the Hamiltonian formulation of Einstein's gravity the canonical variables are the spatial 
metric gy and its conjugate momentum n ij . In our quest for a more general class of theories in 
which the temporal diffeomorphisms are replaced by some other symmetry the canonical variables 
are the same as in Einstein's gravity: = (gij,n ij ). The space I is an n dimensional manifold. 

Let = Ml) be a collection of functionals of the canonical variables. We will consider 
infinitesimal transformations whose generators, which are actually constraints, are ^ fz o^J^, 
where = (a, a 1 ) = (a, a) are the parameters of the transformation, also known as descriptors. 12 
The collection is denoted by a when too many indices would clutter up a formula. Tradition- 
ally, the super-Hamiltonian Jf? and the supermomentum Jtf, called sometimes Hamiltonian and 
momentum constraints, are scalar and covector densities on I of weight 1. a is a scalar, a 1 is 
a vector field. Coordinate independent notations would be cumbersome because of the presence 
of densities of nonzero weight. Any integral in the main text is meant with respect to a fixed 
coordinate volume element on L, and we represent densities t by t = y/gT, where T is a tensor 
independent of the choice of the volume element. On the other hand, in the appendices any integra- 
tion is meant with respect to the natural volume element. Since the pairing of J% with parameter 
a M will appear frequently in our analysis, we will abbreviate it by 

a?-<K ■= Z f cfjf". 

If F is a functional of the canonical variables, its infinitesimal transformation generated by is 
given by the Poisson bracket 

6 s F = {F,a?-Jir M }. 

The parameter functions of and ft 1 have to have some appropriately prescribed boundary values or 
decaying properties in order that the Poisson bracket {a^-J^^-J^}, which is defined by functional 
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derivatives, exist. We shall assume that it exists if both and are in C™, independent of (gy, jf j ). 
The symbol denotes the space of compactly supported smooth functions or vector fields on Z. 

We give our definition of functional differentiability and some other properties of functionals. 
The argument of a general functional is enclosed by brackets, while the variable of an ordinary 
function is in parentheses. 

Definition. Let ^ be a collection of smooth tensor densities on the space Z, V the manifold of 
configurations <f> on which the functionals F : V —> R and H are defined. H[(f>] is a tensor density 
on Z. 

• F is functionally differentiable if for any variation, i.e. for any one-parameter family of 
(<p A )Ae(o,i) e r for which (A, x) i-» <f> A (x) is smooth and <f> A has appropriate boundary values or 
asymptotic properties, the derivative j^F[<f> A ] exists, and there is a collection T[<f> A ] of tensor 
densities on Z such that (A, x) i-» T[(f> A ](x) is smooth, and 

—^r- = Y J T m w A \d A r A , 

aA m L 

where the sum extends over the label m of the members of the collection T and 0, and all the 
tensor indices of T m , which are suppressed, are contracted with the corresponding indices of 
<J> m . The derivative d A $ is simply the partial derivative of the smooth function (A, x) i-> <p A (x) 
with respect to the first variable. T[0] is called the functional derivative of F at 0, and 
denoted by 

d(p(x) 

• H is local in if supp(//[0J - H[<f> 2 ]) c supp^ - <p 2 ) for all <p l 2 e V. 

• H is ultralocal in if it is a function of <f> (but not its derivatives), that is, there is a function 
h such that H[(f>](x) = h{<p{x)) for all x e I. 

• H is concomitant of if = f t H[<f>] for any diffeomorphism / : Z — » Z. 

The definition of "concomitant" gives a precise mathematical meaning to the property which is 
sometimes described as "constructed out of <f> in a form invariant way" or "depends solely on <f>". 
Note the difficulty of giving a sensible definition to this property without assuming any tensorial 
structure on H[(f>] . Hence the definition of "concomitant" includes the assumption of some tensorial 
structure so that the Lie transport of these quantities is defined. 

Unbarred a denotes not only the Oth component of a, but if it is written in place of a (as in 8 a f, 
C^la,^, Q], etc.) then it denotes an a which is given by a = (a, 0). The same convention applies 
to 8$f etc., i.e., if a stands in place of a, it refers to a = (0, a). 

We say that the Poisson algebra of M'^ closes if there are structure functions C M = (C, C) such 
that 

{a»-J%,/F-J%} = C[a,p, g,nVJ%, (8) 

for any ot 1 ,^ 1 e C^°, independent of (gij,n ij ), and a similar equality holds if they are obtained as 
structure functions in some previous Poisson bracket. We indicated that C M are functionals of 
the parameters and the canonical variables. C M [a,/3, g, n] are functions on Z, but that is not the 
reason for their name "structure function." This name refers to their dependence on the canonical 
variable, which prevents the Poisson algebra from having an ordinary Lie algebra structure, in 
which case the term "structure constant" would be appropriate. For further technical assumptions 
see Appendix A. 
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III. Supermomentum 

The supermomentum, is assumed to generate the spatial Lie transport of the canonical vari- 
ables = (gij,n ij ). Let a be a vector field on the n dimensional space L, independent of 0. Our 
assumption on Jtf is that c) 

£ g< f> = 6 3 <f> = {</>, a 1 -m- (9) 
The Poisson bracket of two functionals F and G of is defined as 

kg, (10) 

6gij Stt'J Sgy Sn'J 

Using the Lie derivative of the metric and the momentum, which is a contravariant symmetric 
tensor density of weight 1, we get that J% satisfies 

{ gij , a k -M) = 2V (i of /) = a k d kgij + 2g k{i d j)a k , {n\ a k -M} = d k (a k j j ) - 2n k(i d k a j) . 

If T is a tensor independent of the canonical variables, then the transformation of its covariant 
derivatives is 

S 3 V h . . . V ia T = {V ;i . . . V im T, d-Mh = £*V fl . . . V im T - V,- . . . V im £ s T, 

where we suppressed the indices of T. V is the torsion free covariant derivative compatible with 
the metric gy. Thus for a and $ independent of (gy, 

Lit^gij - -Cp-Cagij = -C&fcgij, and LS = hence 

[Sg, S ] gij = -L m gij = S m gij- (11) 

For any spatial metric we have 

{gy, [dJf-M) = 8 [S 0ij = -[8 S , 8g\gij = -(6 3 60 - 5g8 s )gij 

= -{{ gij ,f3 k -Jm,a l -J%} + {{ gij ,a k -M},/3 l -J^ = {g y M-J? k ,p l -M) 

where in the last equality the Jacobi identity was used. Therefore equation (11) already fixes the 
Poisson algebra of J^: 

[ct-tfuF-'*® = V&fti-Mi (13) 

Note that if a and (3 are independent of the canonical variables, so is [a,fi], hence (9) was indeed 
applicable in the first equality in (12). We quote the following result from [2]: 

Proposition. If is concomitant of the canonical variables (gy,^), and JLsgij = {gij,a k • 3%k), 
SL&ifi = {ifi ,a k ■ ,W k \for any a 1 , ft e C™, independent of (gy, n ij ), then the only possible 3*% are 

^=-2^VsV™ (14) 



c) Strictly speaking, only t-g = J flgy and t-n = f tyrf-* can have well-defined Poisson brackets with the generators 
o^-Jifft, where t is a tensor density. When we write {F[<f>](x), ■ J%^} = G[a, (f>](x), we mean that G^,^] is the 
smooth tensor density for which {t ■ F[<f>], af-Jf^,} = t ■ G[a, 4>] holds for any compactly supported smooth tensor field 
t which is independent of = (g, n). The tensor indices of t, F, and G are suppressed. 
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In section V we will define a dynamics on conformal classes of metrics. To this end, it will 
be useful to realize the spatial diffeomorphism algebra on the reduced phase space of conformal 
classes of metrics. 

Let / : L —> L be a diffeomorphism. A metric gy and g\. are said to be in the same conformal 
class if there is a function <x on I such that g\. = gye a ' . For a generic / the metrics f*gy and 
are not in the same conformal class even if gy and g'~ are. However, it is possible to define the 
transformation of the metric so that it is a map between conformal classes. Define 

^:=^+-VSV,-^=, (15) 

n yfg 

where n = n\. The infinitesimal transformation generated by a' on the metric is 

{gy,a k -m = 2V (iaj) - ^gijV k a k , 

The Poisson algebra of J% is the same as that of given by (13). Let us introduce the conformal 
metric and momentum: 8 

gij '■= g~"gij, & = g" (x' J ~ ^g' Jn j ' ( 16 ) 
The conformal supermomentum (15) is an expression of the conformal variables (16) only: 

J%=-2g ij %7?\ (17) 

where V is the covariant derivative compatible with the metric gy. The transformation of the 
conformal variables generated by a'-Jtf: 

{gy, a k -J%} = 2% aj) - -gyV k a k , {n iJ , a k -Jm = d k (a k n ij ) - 2^d k a j) + -0? V k a k . 

n n 

IV. Super-Hamiltonian 

The question that we turn to now is what are the possible super-Hamiltonians such that the 
Poisson algebra of Jft^ closes in the sense (8), where M'i are given by (14). We shall consider 
super-Hamiltonians which are (i) ultralocal and (ii) at most quadratic in the momentum. We 
also assume that (iii) does not contain a term linear in the momentum n ij , so Jf is the sum 
of a momentum independent potential term and a kinetic term which is a homogeneous quadratic 
expression of n ij . Finally, we assume that (iv) the kinetic term is ultralocal in the metric. 

In [2] condition (i) is the consequence of a kinematical condition: The transformation of the 
spatial metric generated by a- Jt? is required to be ultralocal in a as in general relativity. Hojman, 
Kuchaf, and Teitelboim were interested in the possibility of different Einsteinian geometrody- 
namics, i.e. the existence of such that (J%?, M\) give inequivalent canonical realizations of the 
symmetry algebra of Einstein's gravity. 1 3 If (ii) is also assumed, then Einstein's gravity is the 
only realization of that algebra even for spatial dimensions higher than three. 1 In this case (iii) and 
(iv) are consequences of the symmetry algebra. It was argued in [2, 3] that if the space is three 
dimensional, then even (ii) can be relaxed: The kinetic term is allowed to be a power series of the 
momentum with metric dependent coefficients. In [2] only time reversible geometrodynamics is 
considered, i.e., ^ is taken to be an even function of the momentum, but this condition turns out 
to be redundant. 3 Thus for a three dimensional space (i) is enough to regain Einstein's gravity as 
the only Einsteinian geometrodynamics. 
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Since in our case the symmetry algebra is only partially known, the analysis for such a gen- 
eral class of 3^ as in Hojman, Kuchaf, and Teitelboim's investigations would be much more 
complicated than the exploration of the different realizations of the fixed symmetry algebra of 
Einstein's gravity. Dropping (i) would probably make the analysis of Hojman et al intractable, 
and in our case it would allow for uninteresting modifications of Einstein's gravity. For example 
3^ = (Jif ADM +t ij 3^ i 3^j, J£f), where ^adm is the Hamiltonian constraint of general relativity, and 
is a tensor density of weight -1, concomitant of the canonical variables. Theories with 34? 
satisfying (ii) appear to be the physically most relevant ones since the relationship between the 
momentum and the velocity is linear, so the momentum can be uniquely eliminated in favor of the 
velocity if the kinetic term is a nondegenerate expression of the momentum. Condition (iii) can be 
interpreted as the requirement of time reversibility. 2 For assumption (iv) we do not have any phys- 
ical motivation. It is assumed for the sake of simplicity of our analysis, but we note that Hofava's 
proposal and its later modifications, including the ones in which the detailed balance condition is 
broken, all belong to the class to which the theorem in this section applies. 

The question of different extensions of spatial diffeomorphisms was already raised in [1], where 
Teitelboim notes that we do not know what makes the Poisson algebra of general relativity prefer- 
able over other possibilities, and he mentions ultralocal gravity as an alternative to Einstein's grav- 
ity. As far as we know, all the attempts to find theories which admit an extension of the spatial dif- 
feomorphisms with a further local symmetry, and in which (gy, jfi) are the only canonical variables, 
satisfy condition (i-iv). In [4] the following potentials terms in 34? were tried: V = /uR + v, V = R a ; 
and the potential V = c\R 2 + c 2 RijR ij + cJJ iV l R is also claimed to have been tested (jj, v, a,c\,c 2 ,c^ 
are constants, Ry is the Ricci tensor of the spatial metric gy, R=Ryg ij ). None of them were found 
to result in a closed Poisson algebra of 34?^. In [13] it is conjectured that if the kinetic term in 34? 
is the same as in the Hamiltonian constraint of general relativity, then the Poisson algebra of 34^ 
closes only for Einstein's gravity. Apart from (i-iv), some mathematical properties described in 
Appendix A, and a simplifying condition on the momentum dependence of the structure functions 
that we will describe soon, we do not make further assumptions. In particular, we do not use some 
specific form for the potential, but it is allowed to be a general, possibly nonlocal functional of the 
metric. So our analysis applies to 34? of the form 34? = -^Gyun^n 11 + G[g], where Gyu is ultralocal 
in the spatial metric gy, but G[g] can be a general functional of gy. 

In Einstein's gravity the structure functions C'[a,/3, g] in (8) depend on the metric. In general 
C[or,/3, g, n] can depend both on the metric gy and the momentum n ij . However, we will consider 
only the case when C[a,/3, g, n] is at most linear in the momentum. Recall our notations: This 
means only an assumption on the structure function with which 34? might arise on the right hand 
side of its own Poisson brackets. In the theorem below, for clarity, any possible momentum and 
metric dependence is indicated explicitly. 

Theorem 1. Let the space L be n > 3 dimensional. 34?[g, n] is the super momentum given by 
(14), and 34?[g,n] = G^t^ 1 + G[g], where G m = K(gikgji + gugjk) - Agijgu with constants 
satisfying nX ^ 2k + 0. Assume that for any a,/3 e C~ there is a smooth tensor field t[a,fi, g] such 
that supp t[a,/3, g] c supp a U supp /3, and for any momentum n ij 

C[a,/3,g,n]-34?[g,7r] + C i [a,j3,g,n]-34? i [g,7r] = J ty[a,j3, g]n ij 

i 

with smooth functions C[a,fi, g, 7r] and C l [a,fi, g, n], which are linear in a andfi, and their support 
is within supp a U supp/?. Furthermore, C[a,j3, g, 7i] is also linear in Then there are smooth 
vector fields v[a,fi, g], linear in a and fi, such that supp v[a,fi, g] c supp a U supp/?, and for any 
momentum n ij 

f ty\a,p,gW = f v\a,p,gl^g,nl 
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In other words, under the simplifying assumptions on the structure functions mentioned in 
the above theorem, if the Poisson bracket is a linear expression of the momentum, which is true 
if is of the form as in the above theorem, and it closes, then it is a linear expression of the 
supermomentum itself. So the structure function C[a,/3, g, n] actually vanishes. This statement is 
plausible, but the rigorous proof, which can be found in Appendix E, is somewhat technical, in 
part owing to the relatively weak locality assumptions on C M [a,/3, g, n]. Actually, the statement is 
plausible even if we relax the condition that C[a,fi, g, n] is linear in n ij , especially if we assume that 
C M [a,fi, g, 7i](x) is a polynomial of rt* and its derivatives at x e Z. If the degree of the momentum 
dependence of the structure function accompanying Jf? is as high as in M' itself, the appeal of the 
canonical formalism seems to be lost. We are not going to extend our analysis in this direction. 
The following theorem applies to the case when the Poisson bracket {a ■ J? ', /3 ■ Jf} closes on the 
supermomentum, i.e., C[a,f3, g, n] = 0, which is true if C[a,fi, g, n] is required to be independent 
of the momentum, or if it depends linearly on the momentum, and C M satisfy the conditions of 
the previous theorem. In the next theorem the metric and momentum dependence of Jif^ and the 
structure functions will not be indicated explicitly. The parameters a M ,^ are always assumed to 
be independent of the canonical variables. 

Theorem 2. Assume that the Poisson algebra ofJrf^ = (J$f, J^) closes with the standard Poisson 
bracket (10) on ,3%. J$? M are functions on the n > 2 dimensional space L, concomitant of the 
canonical variables (gy, n ij ), and 

J? = -^pGijunW 1 + G[g], 

£t8v = {gij,af-J®, L^ = WKa l -m for any a\0eC™, 

where G^i is an ultralocal, but G[g] a general (not necessarily local) functional of the metric. 
Under these conditions J4? can always be resettled so that the Poisson algebra of JT^ is 

\a-^^-M'\ = -eg ij (adj]3-{3dja)-J%, 

{ct-Wtf-Wfr = [ctJ] l -J%, 

where 6 = (symmetry algebra of ultralocal gravity) or e = +1 (symmetry algebra of Einstein's 
gravity with Euclidean (e = I) or Lorentzian (e = -I) signature). 

IfGyki is an invertible map between the spaces of rank two symmetric tensors, then G[g] is a 
function of the metric and its first and second derivatives. If e = also holds, then G[g] is actually 
ultralocal in the metric. Ife^O, then Gij k i = K(g ik gji + g a gj k - ^ gygki), where k is constant. 

IfGyki is not invertible, then 6 = 0. 

Proof. The proposition in Section III has already established the only possible form of J#f. We 
divide the analysis of into five steps. See Section II for our notational conventions. 

1. Auxiliary noncanonical variables 

Since Gyki is assumed to be ultralocal in the metric, the Poisson brackets of a-Jti? and ft- can 
produce only linear functionals of the momentum. By the assumption C[a,f3] = 0, these Poisson 
brackets must close on the supermomentum J% with structure functions C'[a,fi] independent of 
the momentum. is assumed to be a scalar density, concomitant of the canonical variables, the 
Poisson bracket of Jt? with fi' ■ J% is the Lie derivative of Jt? with respect to {3, which is a linear 
homogeneous expression of and its first derivatives. C'[a,/3] = 0, and the weight of J$? fixes 
C[a,fi]. What will be important first is that C[a,/3] are all independent of the momentum. Finally, 
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the structure functions C'[a,/3*] = [a,$] 1 , which are independent of the canonical variables, and 
C[a,j$] = are already known from (13). Note that C^{a,p\ is bilinear. 

Let us introduce the auxiliary variables N = N M = (N, N'), which is a collection of a scalar and a 
vector field in CJT, independent of gy and n ij . By the assumption on the momentum dependence of 

<7C , 

\ ... J =2N(x)P ij (x), where Py := -—. 
6tt ij (x) 2 dn lJ 

We introduced the quantity Py for the sake of brevity of later calculations. Note that 

{gy,W-J^} = 2NPy + 2V (i N J> (18) 

By (18) the infinitesimal transformation of gy generated by Jf satisfies 

KSpgv = N{ gij ,/3- J?} = 2N/3Py=p({gy,W-J%} - 2V (i N jl ) . (19) 

In order to get a similar expression for the second variation 5p6 a gij corresponding to two in- 
finitesimal transformations, we have to exchange the order of the transformation and the Poisson- 
bracketing with W-J^. The Jacobi identity tells us that 

{{g^N^},a^} = {{g^a^ 

= [{gy, af-Jty, W-^- [gy, C^[a, N] -J^\ . 

Using that C are independent of the momentum, and a does not depend on the canonical variables, 
we can evaluate the right hand side: 

{{gy,W-^},a-J^} = 2a{Py,W-^} - 2(V (i C /) [a,N] + C[a,N]Pij) • (20) 

From (19) and (20) the second variation of the metric corresponding to two consecutive transfor- 
mations generated by a- and fi-Jf, respectively: 

N6 a 6 l3 gy = {N6 /3 gy,a-jr} 

= 2a/3{Py,N»-J%} - 2j3(V {i C j} [a,N] + C[a,N]Pij + 2aP k(i V j) N k + N k V k (aPy)) , (21) 

where the last two terms came from the evaluation of the Poisson bracket of 2V( ! N 7) = 2g k ^dj)H k + 
N k d k gy with /3-Jf. The commutator of the two transformations: 

N[£ a , S p ]gij = N(S a 5pgy - 5 p 8 a gij) 

= 2Py (a (C[fi, N] + N k dkP) - /3 (C[a, N] + N k d k a)) (22) 
+ 2(aV (i C f) \fi,N] -^V (i C y) [a,N]) . 

Similar calculations lead to the following expression of the commutator of two transformations 
generated by a-Jf and fi'-J^t, respectively: 

mS a , Sfigy = 2NPy/3 k d k a + 2a (Py (C0, N] - /3 k d k N) + g^ (c k 0, N] - £ N k )) . (23) 

2. A note on the tensorial structure of 

is a scalar density of weight 1 by assumption. Multiplying ^ by a power of ^[g, we get an 
of different weight. This freedom in changing the weight of M' is reflected by (22) and (23). 
C\fi, N] = C\fi, N] + C\fi, N]. We have already seen that C\fi, N] = 0. The last term in (23) vanishes 
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since the structure constants in the Poisson algebra of are precisely C l \$, N] = X^N' (see (13)). 
We conclude from (23) that 

C0,N]=0dtN + NC0], (24) 

where C is possibly a functional of the metric, and a linear functional of f3 l . Note that the term 
proportional to Py in (22) completely cancels since C[fi, N] = 0, and C[fi, N] = C[J3, N] = -Ndifi - 
jSC[N] from (24) and the antisymmetry of C in the parameters. The Lie derivative of a scalar 
density of weight w is = 30^) + (w - Since = {J^,0-J^ and 

C0, N] -Jif = C0,N]-J^ + C\fi, N] -{N- JT,/?'- J#?} = -N-£gJ>t? = (pdjN - N(w - l)^')^, 

a scalar density of weight w corresponds to C0] = -{w - l)dfi l . With our choice of weight 1, 

C[ct,/3] = a 1 dfi. (25) 

3. Locality 

Even if we allowed G[g] to be a nonlocal functional of the metric, we can rule out the possibility 
that the structure functions are nonlocal in the parameters. As mentioned in the previous paragraph, 
(22) simplifies to 

m a ,Sp]gy = 2(aV {i C j) lJ3,N] -/JV^C^N]) . (26) 

Also, C[J3, N] = C\fi, N] because C[J3, N] = 0. From (26) we get 

aV (i Cf) \fi,N]\ (supp/3) c = N@ij[a,p]\ (supp/3) c, (27) 

where the superscript c indicates complement. <Py (as C) can depend on the metric. Since are C~ 
can be arbitrary, 

suppV^C^NJI^pp^ c suppN| (supp/3) ,, 

Since the only solution to the Killing equation on an open subset of I is zero for a generic metric, 
we have 

suppC^NJI^upp^ c suppN| (supp/3) ,. 

A set of metrics is not generic if a structure function that is nonzero only on this set could not arise 
in an element of the Poisson algebra since the latter would not be functionally differentiable. By 
property (ii) in Appendix A and lemma 2, there is a finite expansion for any x t supp f3 in terms 
of the derivatives of N so that it gives C,[/3, N] and its first derivatives at x. With this expansion in 
(27), we have 

M 

aV d Z C fl k2 - k '"W {kl V fe . . . V, m) N| (supp/3) c = N0y [<*,/3]| (supp/J) c, 

m 

where the coefficients c k } kl "' km can depend on J3 and the metric, and the number M of the terms 
in the sum may vary with x e L. Since we can always specify N so that precisely one of its 
symmetric covariant derivatives of the highest order appearing on the left hand side is nonzero, 
and all the others, including N itself, are zero at a given point, the only possibility is that all the 
coefficients C k[kl "' km {x) = at x i supp/?. Thus suppC,[/3, N] c supp/3. By property (ii) and 
lemma 2 again, this means that at any point in the entire space I there is a finite combination of the 
derivatives of /3 whose derivatives below a fixed order give C ; and its corresponding derivatives at 
that point. The coefficient functions can depend on the metric and the other parameter N as well. 
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But Ci is antisymmetric in the two parameters, so this expansion can be written in terms of the 
derivatives of the two parameters. So on the entire space I we have 

M 

V, C, \fi, N] = V, (V ft . . . V^/*) V (/i . . . V W N (28) 

m,n=0 

with possibly metric dependent coefficients c kl " Mm ' luJn , which are totally symmetric in k\ . . .k m as 
well as in l x . . . l n , and c kl - km ' h - ln = -C l }~ XnM - km . Note that M can depend on x e I. 

4. Structure functions 

Suppose that at some jc higher than first derivatives arise in (28). Let m(x) be the maximum num- 
ber of ^-indices that occur in (28), and n e N. Fix the following two sets of indices: K := 
{/, k\,..., k m ( X )} and L := {l u . . . /„}. Take an a which is a nonzero constant in a neighborhood of 
x, a /3 for which /3(x) = and dp . . . d k > /3(x) = unless {k'...,k',} = K, and an N such that 

1 m > 1 in 

N(x) = and d? . . . di>N(x) = unless {l' v . . . , l' n ,} = L. What survives from the n-component of 
(26) with this choice is 

= N(x)[6 a , 6 p ]g u (x) = 2a(x) c* 1 "-***' 1 "^*) (d t d k] ... fi{x)) d h . . . ^ N(x), 

where there is no summation over k\ . . . k m ( x ) an d h ■ • ■ h- This can hold for all such choices only 
if n = 0, i.e., highest derivatives of /3 in Ci\J3, N] is multiplied by only undifferentiated N. By the 
antisymmetry of C ( [/0, N] in the parameters, 

C t [J3, N](x) = C k r^\x) (N(x) d h . . . d kmu) p(x) - f3(x) d k ... ^N(x)) + . . . , (29) 

where the ellipsis indicates terms containing derivatives of /3 and N of order lower than m(x). Now 
suppose that m(x) > 2. If we choose an N such that N(x) = 0, and its only derivative that does not 
vanish at x is of order m(x), and f5 is a function such that only dfi{x) ^ for a fixed i, and all its 
other derivatives, including /? itself, are zero at x, then from (26) together with (29) we get 

= N(x)[8a,Sfi\g a (x) = -2C kl - M \x) d^(x) d h ...d km( N(x), 

where there is no summation over k\ . . . k m ( x )- The realization is that it is impossible to satisfy the 
above equality unless m(x) < 1 for all x. So the possible structure functions reduce to 

C i [a,p\ = C{(ad j p-l3d j a) (30) 

with C\ possibly dependent on the metric. By definition, C l [a,/3] and C\ are tensors, but is the latter 
concomitant of the metric? It is hard to imagine that it is not, but let us see a precise argument that 
it indeed has this property. We evaluate the Jacobi identity 

using (25) and (30). We obtain after some algebra that C ij := C' k g jk satisfies 

((y k d k C ij - C ik d k y j - C ki d k i){adfi -fidja)) -J% 

= J d n x {C ij (x), y • m (a(x)djj3(x) - J3(x)d ja (x)) J%(x). (3 1} 

i 

Note that the above equation can be written as f = 0, where £ is a compactly supported 
smooth vector field since a,/3eC™. Integrating by parts, we get that f V (i ^,)^ = 0, which holds 
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for any momentum ', implying that = 0. For a generic metric £ = is the only Killing 
vector field. Furthermore, adfi - fidja can be any vector at a given point. Therefore (31) simply 
means that 

y k d k a - c% r j - c kj d k y = \c\i-m. 

The left hand side is the Lie derivative of a rank two contravariant tensor, so this equation is just 
the mathematical expression of the fact that C ij is a tensor, concomitant of the metric. 

This is not the end of the story. Plugging (30) in (26), and keeping only the terms that are not 
already of the form of an W independent function multiplied by N, we find that 

orC* (djpdft - dfifidp) -pC k (i (dj^dp - d j)a d k N) 

should be proportional to N (with proportionality factor independent of N M ). Let fi be constant in 
a neighborhood of x. Let us choose an N such that N(x) = and all its (first order) derivatives 
vanish at x except for 5,N(x) (for a fixed i), and an a such that d k a{x) is its only nonzero (first 
order) derivative at x. We find that C k = C\8 k (no summation over i). Since C\ was found to be a 
tensor, concomitant of the metric, this implies that 

C[a,p\ = -e[g]g ij (adj/3 - fidja), (32) 

where e[g] is a scalar, concomitant of the metric. 
5. Algebraic properties ofG^i and their implications 

Recall that we chose 3f to be a tensor density of weight 1. Since the momentum is a tensor density 
of weight 1, the coefficient Ggu is a tensor with symmetry properties Gyu = Gjua = Gmj. Gyu is 
concomitant of the metric (this follows from the same property of 34f), and ultralocal in it (by 
assumption). A simple application of Schur's lemma in lemma 7 shows that 

Giju = k (gikgji + gugjk) ~ * gijgki, (33) 

exhausts the class of such tensor densities (k and A are constant). 

Let us evaluate the Poisson brackets in (21), using (25) and (32). Collecting the terms linear in 
the momentum, we get a quantity proportional to 



<>C "" nn V (k N l) - G iMm V n) N k - G mn , (i -V^N*, 



dg, 



ki 



which is zero since the vanishing of this quantity precisely means that the tensor G^u is concomitant 
of the metric and ultralocal in it. The first term is the infinitesimal change of G ijk i if we assume 
that it is a function of the undifferentiated metric, and the infinitesimal transformation is induced 
by the Lie transport of the metric along the vector field N. The Lie derivative of a tensor with the 
symmetry properties of (33) with respect to N gives the other two terms. What remains is 

KSaSpgv = 2Nap[2ptp kl - l-G m ^^G mpq if m iA - 2N/?V (l (eV^) 
\ dg k i yg dgki v? / 

The left hand side is proportional to N, which causes the parenthesis in the second line to be a 
multiple of N, so 

^ NG[g] = V (l (eV / )N) + NFy, (34) 



i 5 r 

—pGiju- — 1 
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where the concrete form of Fy is irrelevant. Recall that we required that the functional derivatives 
are smooth, so Gyki is contracted with a smooth tensor density on the left hand side of (34). Hence 
if Gykiixo) is not an invertible map between the spaces of rank two symmetric tensors at xq, 
then the left hand side cannot span the entire space of SF at x as we vary N. Gyu given by (33) 
is invertible either everywhere {nA ^ 2k ^ 0) or nowhere (nA = 2k or k = 0) on I. In the latter 
case there is no place in the space where the left hand side in (34) generates as N is varied. But 
if e is not identically zero, there are points where the right hand side generates the entire £T since 
the second derivatives of N can be prescribed arbitrarily at one point. The conclusion is that either 
6 = or Gfju is invertible. 

From now on we will consider the case when Gy k i is invertible as a map in the space of rank 
two symmetric tensors. Let G ijkl be the inverse of Gy k i. So it is not obtained by raising the indices 
of G ijkh but by the condition & mn G mnU = l/2(<5j$ + 3$), & u = & m = G ijlk . From (34) we have 

■1=4- f NG[ # ] = G p %(eV N) + NF ij , (35) 

where = & hl F k u F ij = F-*' . Note that (35) implies that G[g] is actually a function of the metric 
and its first and second derivatives (see lemma 8). Using the definition of the functional derivative, 
we can see from (35) that for any variation g A with 6g := d A g A \A=o, 

"tt f NG[gJ = f dgy^- f NG[g]. (36) 
dA J L A=0 J L 8gij J L 

Let the variation g A be given by the diffeomorphisms generated by a vector field We can ex- 
change the order of the integral and the differentiation on the left hand side of (36), plug in (35) 
on the right hand side, and integrate by parts. Considering that N can be any compactly supported 
smooth function, we get that 

6G[g] := d A G[g A ]\A=o = G^eV^V^gy + G^V^V^gy + F*6 gij , (37) 

where we used that the covariant derivatives of G ijkt are zero, and 6G[g] is fixed by the condition 
that it is a scalar density of weight 1, concomitant of the metric. 

Now we argue that the right hand side of (37) is inconsistent with the tensorial structure of 
G[g]. The suspect is the term containing first order derivatives of 8gy. We have 

di(?G[g]) = 2G^eV fe V,V y £ + 2/¥ k V k V^ + 2F%^, (38) 

where we introduced the quantity W k = & kl Vie with the symmetry property W k = W k . Since the 
commutators of covariant derivatives can be expressed in terms of the Riemannian tensor, (38) can 
be written as 

d^G\g\) = 2GWeV ( ,V,V^ + 2#'^%V y) 6 + K"V^ + L%, (39) 

where the concrete form of K ij and V is irrelevant. Since all the symmetric covariant derivatives 
of a vector field £ (including the undifferentiated £ itself) can be prescribed arbitrarily at a point 
(up to a fixed order), and on the left hand side of (39) only £ and its first derivatives arise, (39) 
can hold for any £ only if G^ kt) = and H'W = 0. If e is not identically zero, the first condition 
fixes the value of A in (33): G ijkl = K(g lk g jl + g d g Jk - 2g ij g kl ) with some constant k, and hence 
Gijki = K(g ik gji + gugjk - gygki). Then the second condition yields g^V'e - g i{j V k) e = 0. By a 
contraction we get (n - l)V ! e = 0, so e is a constant. If necessary, rescale J4? by -j= to get e = +1 
in the Poisson algebra of J^ M . 
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Finally, let us assume that 6 = 0, and Gyu is still invertible. For brevity, we introduce the tensor 
density t ij [g,a] = /-(cr-G[g]). The assumption e = implies that 




The momentum n ij can be any compactly supported symmetric tensor. The same applies to GyuTf 3 
by the invertibility of Gyki. Thus (40) implies that t ij [a,g] = aT i] \g\. By lemma 8, G[g](x) is a 
function of gy(x), or in other words, G[g] is ultralocal in the metric. □ 

Now we can decide if Hofava's Lagrangian, or its generalization to any (even nonlocal) poten- 
tial term, is a viable modification of Einstein's gravity if we demand that such a theory should 
give rise to a local Hamiltonian constraint J4f which - under our conditions on the structure 
functions - forms a closed Poisson algebra with the momentum constraint. In general relativ- 
ity G m = K(g ik gji + gngjk - gijgkd- In the class of theories we are considering the coefficient 
in the kinetic term is G^ t = K(gikgji + gugjk) - ^gijgki with some constants k and A. Since the 
goal is such a theory that interpolates between Einstein's gravity and a yet unknown UV theory, 
A arbitrarily approaches ^ k, which means that G^ has to be invertible in some regime of the 
parameters characterizing the interpolation. The following result applies to this case: 

Corollary. Let the space L be n > 3 dimensional, and Gyu is an ultralocal concomitant of the 
metric, invertible as a map between spaces of rank two symmetric tensors. Assume that the Poisson 
algebra of the supermomentum J^i and 

= —G^,n ij 7T kl + G[g] 

closes with the standard Poisson bracket (10) so that the super-Hamiltonian arises in its own 
Poisson brackets with structure function C[a,j3, g, n] which is either independent ofn'i, or it de- 
pends on it at most linearly, and in this case the partial locality condition supp C^[a,fi, g,n] c 
suppa U supp/3 holds for all the structure functions. Then G[g] = ^JgQxR + v), where p and v 
are constant, and R is the Ricci scalar ofgy. If p. £ 0, then G ijk i = k (gikgji + gugjk - ^ gijgki), where 
k ± is constant. 

Proof. According to the theorem, if e = and Gyu is invertible, then the potential term in M' is 
ultralocal in the metric. So G[g] is a scalar density of weight 1, ultralocal in the metric, hence it 
is p Vg. Indeed, ar%^=G[g] = 2(^-^G[g])v,^ for any a\ which implies that V,--^ G[g] = 0. If 
6^0, then by rescaling J$? by a constant, the Poisson algebra, if it closes, can be brought into the 
symmetry algebra of Einstein's gravity. As shown in [1], the only realization of this algebra with 
Jtf given by (14) and an J$? which is a quadratic function of the momentum is Einstein's gravity. 
6 ^ if and only if p ^ 0, and the theorem gives the form of Gp/ in this case. □ 

V. Conformal Lifshitz gravity 

The main motivation of the previous chapter was to see if there are relatively simple theories that 
can be considered as alternatives to Einstein's gravity. Our approach was rather conservative. The 
theory is a geometrodynamics, meaning that in the Hamiltonian formalism it describes the time 
evolution of spatial geometries, and the Hamiltonian is given by a homogeneous linear expression 
of local constraints whose Poisson algebra is closed. The local constraints follow from the Hamil- 
tonian equations. We were interested in an extension of the algebra of the spatial diffeomorphisms 
by an additional local symmetry mostly because the possibility of a continuous deformation of 
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the symmetry algebra of general relativity is tantalizing. In our view, Hofava's attempt to recover 
Einstein's gravity as a limit of some flow of theories can be successful only if such deformations 
exist. Furthermore, with the additional constraint function and the nontrivial transformations it 
generates on the constraint manifold, the naive counting of the degrees of freedom at a spacetime 
point yields the same number as in Einstein's gravity. This is reassuring if the degrees of freedom 
in Einstein's gravity are to be matched with those of the modified theory. On the other hand, as 
already mentioned in the Introduction and Section II, the symmetries associated with the Hamil- 
tonian constraint make it possible in Einstein's gravity to eliminate the metric mode with negative 
kinetic energy. If this is considered to be the key property, one can take a less conservative stand- 
point. The goal can be a theory such that this mode is absent. The number of degrees of freedom 
per spacetime point is then the same as in general relativity, the dynamics is the time evolution of 
spatial geometries, but the Hamiltonian is not necessarily a combination of local constraints. 

As a starting point, we investigate the Lagrangian (7) suggested by Hofava as a candidate for 
the UV fixed point of his theory under the anisotropic scale transformation (in n = 3 dimensional 
space). The properties of this model lead us to the consideration of theories in which the only 
canonical variables are the conformal metric and momentum (16), and their transformations under 
spatial diffeomorphisms are generated by the conformal supermomentum (17). 

In (7) the kinematically possible configurations of N cannot be restricted to spatial constants 
since the scale transformation of N proposed by Hofava is space dependent. The variation of 
(7) with respect to N gives a local constraint which does not form a closed Poisson algebra with 
the conformal momentum constraint. The simplest way of getting a closed Poisson algebra of 
constraints is to demand that N is a (spatial) constant, so the only local constraint is the conformal 
supermomentum. We show that this restriction is not an obstacle to the realization of the local 
Weyl symmetry. We will modify (7) so that the action in the Hamiltonian formalism depends only 
on the conformal variables. Local scale transformations act on these variables trivially, so they 
leave the action invariant. 

The action (7) can be got from the following action if we eliminate n ij , using the field equation 
obtained by the variation with respect to n ij , with the assumption that no spatial boundary term 
arises: 

Suv = / (*% - - N^ uv ) , (41) 

Rxl 

where we used the conformal variables (16). For simplicity we set k 2 = 2. Jtf is the conformal 
supermomentum (17), and 

J«6v = -L(»»»j-i^) + ^C*C,. (42) 

The first term in (41) is a nondegenerate bilinear expression of gy and ffi, so the canonical variables 
are the conformal metric and momentum. This term has the standard form J L n ij gij for which the 
Poisson brackets are defined by (10), but here (gij,n ij ) are subject to the second class constraints 
g = 1 and n = 0, so the bracket is the corresponding Dirac bracket. Since g and n have vanishing 
Poisson brackets with (gy, n ij ), the Dirac bracket is the standard Poisson bracket written in terms of 
the conformal variables and (g, n), without the term containing functional derivatives with respect 
to the latter variables. 

From the perspective of Dirac's approach to constrained systems (see for example [8]) it would 
be more natural to include n, the generator of the local scale transformations of the spatial metric, 
in the Hamiltonian as one of the constraints: 

Suv = / (fgij ~ WJZ-Att) -fH, (43) 

Rxl R 
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where A is a scalar function. Since n is spatial scalar density, the Poisson algebra of (n, Jiff) closes. 
We also included a term H. If the Poisson brackets of H with the constraint functions J£f and n are 
their homogeneous combinations, there is no need for imposing further local constraints in order 
to guarantee that the time evolution determined by (43) preserves n = and 3% = 0. Tentatively, 
we can write H = N-J%y. All the symmetry generators are already present in the Hamiltonian of 

(43) , so we do not lose anything, but gain a closing Poisson algebra of the constraints (n, J^) if we 
modify the action (43) so that is not a local constraint. First we replace the space dependent 
N by a constant v, otherwise the variation of (43) with respect to N would give J^uv = 0. But 
this is not enough in order to avoid ^fj V = as the consequence of the field equations. Using the 
scaling properties of the Cotton tensor (see below equation (7)) algebra, we can see that 

{a-J^ uy Mx)} = ~a(x)jr uy (x), (44) 

which means that n = forces us to impose the Hamiltonian constraint Jtfjy = as well, otherwise 
n would not be preserved by the time evolution. (A negative constant can be added to J^uv in order 
to avoid the constraint n ij = 0.) We can exclude the possibility that the Poisson bracket of a-J^jw 
and fi • closes on the supermomentum. Indeed, if that was the case, then by theorem 2, this 
Poisson bracket should be zero for any a and /3 since Gy k i is not invertible. On the other hand, 

(44) implies that (J#uy := ^uv + ^ would form a Poisson algebra with the same structure 
functions as those of the algebra of (J^jy, Jtf). But the coefficient in the kinetic term of J^jy is 
invertible, so by the theorem again, the potential term in J^jy should be an ultralocal functional of 
gy, and 0->Cij is not of this form. This argument applies to any Weyl covariant nontrivial potential. 

If our goal is such a theory in which the temporal diffeomorphism symmetry is replaced by 
a local scale transformation of the spatial metric, then there is a simple way to save (41). Since 
^uv = is the field equation that we get by varying (43) with respect to the local scale factor of 
the spatial metric, we will eliminate this degree of freedom by imposing the constraint In g = 0. 
The action is 

S uv = / (n ij gij - N'JZ -An-ylng)- f v ^ uv , (45) 

Rxl Rxl 

where v is a (spatial) constant. The field equations of (45) after the elimination of A and y are 
equivalent to the Hamiltonian equations of 

tf uv = ~ 2 f L + v f L (*% + T^q) (46) 

with the standard brackets (10) replaced by the Dirac brackets for the constraints n = 0, Ing = 0. 
V is the covariant derivative compatible with gy, C ij is the Cotton tensor of gy, and the indices are 
lowered by gy. The scale transformation acts trivially on the conformal variables (16). Finally we 
note that naive counting of the degrees of freedom of (46) at a spacetime point gives the same result 
as in Einstein's gravity. The number of the canonical pairs is less by one, but the same applies to 
the number of constraint functions and gauge fixing conditions. 

VI. Discussion 

In Einstein's gravity, the ultralocal metric dependence of the kinetic term in Jf is the consequence 
of the Poisson algebra of M 3 ^} Allowing the kinetic term to depend on the derivatives of the metric 
might open the way to new extensions of the spatial diffeomorphisms. Probably, the relaxation of 
time reversibility is the most interesting generalization. This would mean a term in 3tf linear in the 
momentum. 2 As we mentioned earlier, there is no irreversible Einsteinian geometrodynamics. 3 It 
is an intriguing question if irreversibility allows for a symmetry algebra other than that of general 
relativity or ultralocal gravity. 
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Appendix A: Notations, assumptions, and their immediate consequences 

The space L is an n dimensional smooth Riemannian manifold with a metric g on it. It is not 
required to be compact, and it can have a boundary. Let CI c I be a bounded domain (i.e. a 
connected open set whose closure is compact). The boundary dCl of CI is always assumed to be a 
smooth manifold. C°°(n) is the space of smooth functions on CI, and its subspace C™(Cl) consists 
of functions whose support is within CI. The space of smooth vector, m-index, or symmetric two- 
index tensor fields is denoted by C°°(C1, T), C°°(0, ® m T), or C°°(0, TVT). We introduce the spaces 
C™(C1, T), ... similarly to the scalar functions. C°°(L) is the space of the smooth functions on L, 
and CT(I) is the subspace of functions whose support is in the interior of I. 

The p-integrable functions, vector fields, . . . live in L P (C1), L P (C1, T), ... For p = 2 they are 
Hilbert spaces over R. H e (Cl), H C {C1, T), . . . are the Sobolev spaces, and Ho(Cl), H C {C1, T), . . . are 
the completion of C™(C1), C™(C1, T), ... in them. These are also Hilbert spaces, but we will use 
only the norm, and never the scalar product on them. Thus (, > always refers to the L 2 scalar 
product, defined as (s, t) := f n s ll - lm t il .j m for t, s e L 2 (C1, ® m T) even if the arguments happen to be 
in H c c L 2 . Recall the definition of the Sobolev spaces H e (Cl). (The other spaces H e (Cl, T), . . . are 
defined in the same way.) H°(C1) = L 2 (C1), and the elements of H e {Cl) {£ > 1) are the functions in 
H e ~ l (Cl) whose weak derivative of order t exists, and it is square integrable, i.e., it is in L 2 (C1, ®(T). 
The standard norm on H e (Cl) is defined by ||/||^ (Q) := ||/|&,_ 1(Q) + II t^Sf\\ 2 L 2 (n ^ Ty 

£ 

We will also use the Banach spaces C r (CX) of r times continuously differentiable functions. For 
r = this is the space of the continuous functions, which is usually denoted by C(C1). As usual, 
the symbol C(C1, T), ... stands for the vector and tensor fields. The norms on these spaces are 
ll/llc(0) = max max | f^.^f(x)\, etc., where \t(x)\ can be \t(x)\ 2 = t il .. Jm (x)f 1 - im (x). 

(K.v<r xe ci v. 

s 

The valence of the elements of all these spaces must be fixed, otherwise addition on them would 
not make any sense, but a field whose valence is different from our convention can be considered 
as an element of our space by raising or lowering the appropriate indices. We can define a tensor 
field p by yfgp = n. We will use p instead of the momentum 7r so that we have to deal with only 
tensors, and not with tensor densities as well. 

A kinematically allowed configuration of the canonical variables is a configuration for which 
the (closed) Poisson algebra generated by { a M ■ Jif^ | a M e C™(L) } is defined. (a p are independent 
of the canonical variables.) Apart from smoothness requirements, these configurations are subject 
to some boundary conditions, or they have specific decaying properties if I has a boundary, or it 
is not compact. We do not need a detailed definition of their manifold, but we assume that locally 
the metric and the momentum can be anything, and the manifold of the metric configurations has 
a kind of topological property that makes it possible to infer local metric dependence from the 
functional derivatives of certain functionals: 

i. Let F = Q x V be the manifold of the kinematically allowed configurations of the canonical 
variables, where Q is the space of the kinematically allowed metrics, V is that of the mo- 
menta. C"(L, TMT) c V is a vector space. Let gj and g 2 be two metrics on I, gj 6 Q, and 
g\- g 2 ^ C™(L). Then we assume that g 2 e Q. Furthermore, any g x , g 2 e V are connected 
by a variation (0^)^ e (o,i) e T which is constant on I \ supp(^j - g 2 ). 
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When we say momentum, we always mean a smooth tensor field p which is kinematically allowed, 
pep. Sometimes the word "momentum" refers to the corresponding tensor density n = ^fgp. 
(The momentum conjugate to the metric is actually a tensor density.) The momentum p is said to 
be transverse (on II) if V t p ij = (on CI). 

We will impose some continuity conditions on the quantities arising in the Poisson algebra. 
Let CI c I be open. S>(C1) is the vector space of the compactly supported smooth functions on 
CI with the topology in which a sequence f k goes to zero if there is a compact set K c CI such 
that supp f k c K for all k, and f k and all its derivatives uniformly tend to zero. A sequence f k is 
said to converge to / in @>(C1), if f k - f goes to zero in 3)(Cl). The spaces S>(C1, T), ... and the 
convergences in them are defined similarly. The conditions: 

ii. Let CI c I be open. If (a,p) i-» C[a,/3, g,p]\ a is a map C™(£>) x Cf(L, TwT) —> C™(C1) 
for all g € Q and fi e C£°(L), then it is required to be continuous in the ^-topology. 

iii. Let F : Q x V —* R be an element of the Poisson algebra. Then g \-> F[g,p] is required to 
be continuous in the ^-topology for any pep. 

iv. If g is analytic in local coordinates (CI, cp), so is the potential function G[g] of Jf. 

The space L is defined to be a smooth, but not an analytic manifold. That is why analiticity of func- 
tions and tensors is a coordinate dependent property, as in the last condition. These assumptions 
are not restrictive at all for practical purposes. If G[g](x) is a polynomial of the metric, its inverse, 
and its derivatives at x e L up to an arbitrarily high fixed order with coefficients independent of 
x, then a-G[g] satisfies conditions (iii) and (iv) for any a e C™(L) that does not depend on the 
canonical variables. For theorem 2 we need only the first two properties. The other two are used 
only for theorem 1. Condition (iii) makes it possible to extend certain results obtained for locally 
analytic metrics to any kinematically allowed metric since we have 

Lemma 1. Let Cl$ be a bounded domain for which there is a coordinate chart (CI, ip) so that 
Cl c CI. For any f e C°°(L) there is a sequence f k which is analytic on C1 Q in the given local 
coordinates and tends to f in S>(L) as k —* oo. 

Proof. We can assume that CI is bounded. There is a compactly supported smooth F : R n —> R 
which is equal to foy~ { on tp(Cl). The norm defined by ||F|| 2 = J^„ d"x(l + \x\ 2 ) e \F(x)\ 2 is equivalent 
to the standard norm || • \\ H ((w), where F is the Fourier transform of F on the Euclidean space R". 
Since \\F\y m < oo for any t > 0, lim*^ [ A>k d n x(l + \x\ 2 ) l \F(x)\ 2 = 0. Let F k = F(x) if 
\x\ < k, and F k (x) = otherwise. F k , the inverse Fourier transform of F k , is real analytic for any 
k. There is a smooth function p : CI — » R which is 1 on Cl and supp p c CI. Define f k on CI 
by fk = p(F k o ip) + (1 - p)f, and let f k = f otherwise. By the Sobolev embedding theorem, if 
2(1 - r) > n, then H e (Cl) <^-> C(C1), and this embedding is continuous. Thus the sequence f k has 
the desired property since all its derivatives uniformly converge. □ 

The lemma we just proved naturally extends to vector and tensor fields, including nondegener- 
ate tensor fields like metrics. One example why (ii) can be useful: 

Lemma 2. Let CI c L be open, H : 2l(d) — > 2l(d) a continuous linear map, and supp H[a] c 
supp a for any a e C^(C1). Then at any point x e CI there are smooth functions C J ' lll2 ~ lm so that 
H[a](x) and its derivatives djH[a](x) (where j is a multiindex) can be written for all |j| < /< oo as 

M J (x) 

d i H[a](x)= J] d J (C J ™" J ~d h d h ...d im a)(.x), 

m=0 

where M J (x) < oo for all x and J. 
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Proof. Define 

(T x ,a):=H[a](x), a e C^O). 

By the assumption on /, T x is a distribution. Furthermore, if x t supp a then x g supp //[or], so 
(T x ,a)=0. Thus supp T x c {x}. Any distribution whose support is (at most) a point is a (finite) 
combination of the Dirac-delta and its derivatives at that point: 

M°(x) 

(T x ,a) = J] (-l) m C hh -Hx)(d h d h ■ ..d in 6 x ,a) 

(47) 

M°(x) 

= J] C hh ~Hx) d h d i2 ■ ■ ■ d im a{x), M\x) <oo, ae C?(C1). 

m=0 

Note that we have not proved that M°(x) is bounded in some neighborhood of any x, so we have 
not shown that any derivative of H[a] can be got by differentiating the formula for j = 0. But 
we do not need this. As a &{C1) —* &{C1) map, any derivative of H[a] has the same properties 
as H[a] itself, so it can be written as (47): djH[a](x) = Y^Iq Cj"' 2 -'"' (x) d h d h . ..d im a(x), where 
Mi(x) < oo. djH[a](x) contains only a finite number of derivatives of a. In a neighborhood of 
x let a 6 C™(C1) be the Taylor series of a about x truncated at M J (x), which is the order of the 
highest derivative of a arising in djH[a](x), |j| < J. Clearly, d^H[a](x) = d^H[a](x) if |j| < J. But 
the sum extends over m < M J (x) in a neighborhood of x if it is evaluated on a. The coefficients 
Qi\-im are smooth since H[a] is smooth for any a, and djH[a](x) is simply the derivative of the 
formula for H[a], which is a finite combination of the derivatives of a with smooth coefficients in 
a neighborhood of x. □ 



Appendix B: No hidden constraints from spatial diffeomorphisms 

Let S g be a subspace of compactly supported smooth symmetric tensor fields on the space I. 
Subscript g indicates that the S g might depend on the metric. Suppose that S g is orthogonal to all 
the transverse momenta. In other words, the constraints f L typ^ (t e S g ) are consequences of the 
momentum constraint: = implies that J L tijp' j = for all t 6 S g . (The tensor p is related to 
the canonical momentum by n = yfg p.) Is it possible to express the latter constraints in terms of 
the supermomentum? Is there a smooth vector field v for any t e S g such that f L tijp ij = f L VjVjp ij , 
and the support of v is within that of t? Let the formal adjoint of the differential operator Vip ij be 
K, which is proportional to the Killing operator defined on compactly supported vector fields. It 
would be enough to show that S g c ran K. This is not obvious at all as the relationship ker K + = 
(ran K) 1 - might suggest. One problem is that we should define the Hilbert space in which the 
orthogonal complement, the adjoint, and later the closure are taken. Furthermore, K + is not a 
formal adjoint, but the actual one, and it is likely to be defined on a space bigger than that of the 
smooth vector fields, so it might be an extension of V,//^. Thus (ker K+) 1 - might be smaller than 
S g . Finally, what we have is (ker = ran K, and it is not obvious why the smooth tensor fields 
in the closure ran K should be in fact in ran K. These are the details that are worked out in the 
lemma in this appendix. But before the lemma, we need to ask our question in a more precise way. 

Later it will suffice if certain constraints that follow from the momentum constraint can be 
expressed in terms of the supermomentum only locally, so first we reformulate our question ac- 
cordingly. Let CI c I be a bounded domain. Define the linear functional L : C°°(n, T V T) —* R 
by L[p] := f a tijp ij , where t e C^(C1, T\/T). Let D : C°°(C1, TvT) C°°(D., T) be the supermo- 
mentum considered as a linear functional of p, that is, D-*[p] = Vip'-*. If ker D c kerL, is there a 
v e C™(C1, T) such that L[p] = f n ViVjp^l This is the relevant question in the context of the paper, 



22 



but as it frequently happens, it is easier to give an affirmative answer to it if we rephrase it on more 
convenient spaces, which are larger than those of the smooth functions. We will impose on L a 
kind of continuity condition with respect to the operator D. The expression L[p] := f a typ'-* defines 
a continuous functional on L 2 (C1, TvT). The operator D, considered as a map from L 2 (C1,T\/T) 
into L 2 (C1, T), is defined on a dense subspace, and its closure D exists since the domain of D + 
is also dense in L 2 (C1, T). For example, dom D + surely contains C™(C1, T), on which D + acts as 
£\->-l/2 £^g. Does ker D c ker L imply that L[p] = f L ViVjp ij with some v 6 C™(C1, T)l 

The additional continuity property of L means that we not only require that L is zero whenever 
D is zero, but L also vanishes on tensors which are the limits of sequences of smooth tensor fields 
along which D goes to zero. Note that it is important that v whose existence our question addresses 
is required to be supported within supp t. The following statement will suffice for our purposes. 

Lemma 3. Let £1 be a sufficiently small bounded star domain with smooth boundary dCl. The 
map L : L 2 (C1, TvT) —> Ris defined by L[p] = f n typ^, where t e C™(C1, TvT), and the operator 
D : L 2 (Q, TvT) >-> L 2 (C1, T) is given by domZ? = C°°(0, TvT), D\p\ = V iP \ If 'ker D c kerL, 
then there exists a v e C™(C1, T) such that L[p] = f n ViVjp ij for any p e C°°(C1, TvT). 

Proof Let K : L 2 (C1, T) >-» L 2 (C1, TwT), dom K = C™(C1, T), K[£] = -l/2£gg. We have K + D 
D since K + is closed, and K + \ c ^ { Q_ TyT) = D. (A d B means that dom A d dom B and A| domB = B.) 
Actually, K + = D. The argument is borrowed from the theory of Sobolev spaces. The existence 
of the derivative w-* = V,/^ in the weak sense supplemented with the square integrability of p and 
w is equivalent to p e dom K + . Along the lines of the argument that shows that C°°(n) is dense 
in the Sobolev spaces H e (Cl) (this is where the star-convexity of CI comes into play), here we get 
that the elements of dom K + are the limits of convergent sequences p k e C°°(C1, TvT) for which 
Vip'l is also convergent. This precisely means that dom K + c dom D. 

Since k&x K + = (mnK)- 1 , where J_ denotes the orthogonal complement, (ker K + Y = mnK, 
and therefore (ker D) 1 - = mnK. By the assumption on L, the tensor field t is orthogonal to 
kerl), hence t e ran K. What we have to prove is that in fact t e ran K, or in other words, that 
ranK r\_C™(fl, TvT) = ran K. 

Let & e CT(I1, r) be a sequence such that K[g k ] is convergent in L 2 (D.,T\/T). First we want 
to deduce the convergence of & from that of £g k g. The Riemannian version of Korn's inequality 14 
states that there is a constant C such that 

IIV^I 2 < C(||#|| 2 + HX^II 2 ). (48) 

for any £ e C°°(D., T) (and hence for any £ e H l (£l, T)). (The condition on dD. can be weakened. 
See [14] for the details.) The symbol V£ denotes the rank two tensor of the co variant derivative of 
By Friedrichs' inequality, 

ll^llo < (diamo 0)||V #||o (49) 

for any for any | e C™(C1, T) (and hence for any £ e Hl(Cl, T)). The subscripts in (49) indi- 
cate that the norms and the derivatives are calculated on the flat Euclidean background. Note the 
boundary condition on £ for Friedrichs' inequality. The coefficient diam CI is the diameter of 12, 
that is, the longest distance between two points of the boundary. The boundedness of the metric 
and its derivatives on O allows for the following straightforward, even if perhaps not the most 
efficient generalization of Friedrichs' inequality to a Riemannian manifold: 

||^|| 2 <c(diama) 2 (||V|| 2 + ||^|| 2 ), (50) 

where c is a constant. This and the other constant C in (48) can depend on g and CI. The derivation 
of (50) from (49) also shows that c can be chosen so that it depends on the metric g only through 
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a continuous function of the maximum of the Chnstoffel symbols on CI (maxQ r'.£g jm g* cn r^ mn ), 
max xeci g( x ), and 1 / niin xe Q g(x), where the extrema of the matrices in a given coordinate system 
are the minimum and the maximum of their eigenvalues on CI. Let Cl be a bounded domain with 
smooth boundary, and assume that all the domains CI we are considering here are within Cl . If 
£ e C™(C1), then the extension ft is in C™(C1 ), where ft| a = ft and fth \n = 0. This means that 
similarly to c, the constant C can also be chosen to be Q-independent so that (48) and (50) hold for 
any CI c Cl . Thus the combination of these inequalities results in an upper bound on the Sobolev 
norm ||£||//i(n,r) of £ e C™(C1, T), provided that CI is sufficiently small: 

lliw) := ^ l|2 + H^H 2 < rllX^II 2 , 

where the constant y depends on the metric and the domain. Since £,^g converges in L 2 (C1, TvT), 
so does ft in H l (Cl, T). 

Next, we will show that £ := lim^ft e C°°(Q, T). Since the limit is taken in H l {Cl, T), g 
is also in H l (Cl, T). In particular, the first weak derivatives of £ exist, they are the limits of the 
classical derivatives of ft in L 2 (C1, T <g> T), so we have V,-^ + \^-ft = -ty, where the derivatives of £, 
on the left hand side are defined weakly. Thus 

<V ; .v ; Uv> = </./,;) -(V,^ 7 .^) forany a e C?(C1,T ® T). (51) 

The smoothness of | can be shown by induction. Let be the Riemannian curvature tensor. 
Following the steps by which V ; -V ; -ft = -Rjk/^i is derived for a Killing vector field, we can prove 
that the second weak derivative of £ exists. We start with 

<(V,-V 7 - - VjVds ijk , ft) = ft) for any a e C?(Cl, T ®T ®T), 

where (, > is the L 2 scalar product. Using (51) for replacing the second term on the left hand side, 
we get 

(V.-V//* + s jki ), ft) = -</* V, ft) + (V,/*, ^>. 

Now, we can write down this equation for other two tensors, a' and a", whose components are 
permutations of those of s, namely, s' ijk = Sju and s"- k = Skij. Adding the s'-equation to the s- 
equation, and subtracting the s"-equation, we get 

(V,- V*, ft) = R k ^ - V itjk - V jtki + V^>, (52) 

which means that the second weak derivative VV£ e Lr{Cl, T ®T ®T) exists: 

V;V,-ft = R kij % - - Vjt ki + V k tij, 

so £ € H 2 (C1, T). Note that the product of a weakly differentiable and a smooth function is weakly 
differentiable. Let us replace s ijh by Vis ijkl in (52). We can recast the derivatives of s ijkl on the other 
argument of the scalar product on the right hand side because the lower order weak derivatives 
of £ have already been proved to exist. Hence we conclude that the third weak derivative also 
exists, moreover, £ e H 3 (C1, T). Moving on to the higher derivatives in this manner, we can 
see that ^ e C\fy\H e {Cl, T). In order to conclude smoothness from this, we apply the Sobolev 
embedding theorem, which states that if 2(£ - r) > n, where n is the spatial dimension, then 
H e (Cl, T) <^-> C(C1, T), and this embedding is compact, in particular, continuous. 

What remained is the proof that not only | e C°°(C1, T), but also | e C™(C1, T). At this 
point we could refer to the nonintegrability of the Killing equation for a generic metric, but Korn's 
inequality will be enough. If n = 1, then H l {Cl,T) <^-> C(C1,T), implying that the convergence 
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of % k is uniform, so it is also pointwise, therefore g\ m = 0. For n > 1 the condition $\ dn = 
does not follow from the convergence of g k in H l (Q, T). But recall that the convergence of g k and 
K[g k ] means that £ e dom K. Recall that K + = D, so we have K = K ++ = D + (= D + ). Hence 
/? i — * (£, -D[p]) = fn£iVjP ij = § m ^(inj)p ij - f n V^j)p i] is continuous. The vector n is the outward 
normal to dCl. Since p\ga can be any smooth function so that \\p\\ is arbitrarily small, we have 
^(iHj) = on dCl, which implies that %\ aa = 0. 

Since dD. is smooth, there are coordinates in a neighborhood U of x such that dCl is an n - 1 
dimensional plane. Let to be a like a mercury droplet sitting on a horizontal pane of glass. That is, 
u> c [/ n Q. is a domain with a smooth boundary, and there is a ball 5 of nonzero radius such that 
do) n dCl n 5 ± 0. For smooth vector fields on that vanish on do> n 5, Korn's inequality holds 
without the term proportional to ||£|| 2 in (48). 14 Since supp t c O, if co is sufficiently small, t\ w = 0. 
Thus V|L = 0, which - together with g\ m = - implies that £L = 0, therefore £ e C~(a, T). □ 

Appendix C: Linear elasticity and the momentum constraint 

Let us fix the metric, and say that we want to solve the momentum constraint V,/? 1 ^ = for p in 
some domain 12 of the n dimensional space. Is there some simple expression of p that can be 
freely prescribed apart from certain regularity conditions which come from the requirement that 
p e H l (Q.,TvT) or p e C°°(D., TvT)7 The simplest guess, the trace of p is actually a promising 
candidate. The resulting equation is the equilibrium equation in linear elasticity: V,-cH' + f = 0. 
The traceless part of p plays the role of the unknown stress tensor <x, and the force density / 
corresponds to 1 /n Vp, where p = p\ is the prescribed trace. In the Euclidean space the convolution 
of / with Thomson's solution gives the deformation u in W\ The Thomson kernel can be used to 
prove that there exists a smooth solution p if p is smooth. In a curved space we do not have the 
luxury of having such a fundamental solution at hand. But Korn's inequality, which was initially 
motivated by linear elasticity, is helpful even in this case. The original inequality underwent several 
improvements: not only flatness of the metric has been relaxed in [14], but also the strain - the 
Lie derivative of the metric along the deformation - has been replaced by its traceless part in [15]. 
This is the relevant inequality in this section. 

As usual, it is easier to prove the existence of a traceless solution S[f] to V,-S^[/] + f J ' = 
if the derivative ViS ij [f] is defined only weakly. Here the symbol S[f] just emphasizes that we 
are looking for an assignment of solution to the source /. This assignment preferably has some 
continuity property. The regularity on / can also be lowered if the differential operator is only 
weakly defined, and ultimately we are seeking an everywhere defined continuous solution assign- 
ment S : L 2 (C1, T) —* L 2 (Q, TvT). If there are nontrivial (weak) solutions p to the homogeneous 
equation V t p ij = 0, then this assignment is not unique. The existence of such nontrivial solutions in 
a bounded domain is essential in linear elasticity. This is the freedom that makes various boundary 
conditions imposed on the stress and the deformation admissible. The most natural way to resolve 
this ambiguity is to project out the solutions to the homogeneous equation. We demand that S[f] 
be orthogonal to all the weak solutions to Vjp ij = 0. 

This choice is natural in linear elasticity as well since it corresponds to a simple boundary 
condition on the deformation u, as we will see soon. However, the elastic material whose equilib- 
rium equations are analogous to the momentum constraint is rather exotic. As noted in [15], the 
corresponding equations of linear elasticity can be derived by varying the energy E with respect 
to u, where the energy is defined by AE = ||Xjtg|| 2 - A(f,u) in which £,ng := £,ug - 2/ngVjU l 
is the conformal Killing operator. In the Euclidean space this energy is invariant not only under 
isometries, but also similarities. That indicates that the bulk modulus of the material is zero. As far 
as we know, there is no elastic material that is infinitely compressible, but resists shearing. Nev- 
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ertheless, we can imagine that this fictitious material is glued to the boundary dCl of the bounded 
domain 12. Then the equilibrium equation reads as V ! (Xj?g) 1 -' + f J = with the boundary condition 
m|«9q = 0- Thus p = £,ag is in the image of the conformal Killing operator defined on vector fields 
vanishing on the boundary, hence it is orthogonal to the kernel of the adjoint of this operator, which 
is nothing but the weakly defined differential operator Vip lj with no boundary condition on p. So 
if we accept the idea of our exotic elastic material, we have no physical ground for doubts as to 
the existence of the map S with the specification that ran S is orthogonal to the (weak) solutions 
to the homogeneous equation. The upcoming lemma gives a proof for this expectation. 

Finally, we note that it is advantageous to introduce the deformation u even in the investigation 
of the momentum constraint, where it is only an auxiliary variable. The resulting second order 
equation V^-Lag)^ + f J ' = for u is elliptic, so the arsenal developed for elliptic equations can be 
deployed. 16 Alternatively, one can use a generalization of the Lax-Milgram theorem after proving 
that its conditions hold, for which Korn's inequality can be used. 17 We will not follow these paths 
here, since we do not need so strong results as the ones that can be obtained by these techniques. 
Our existence theorem is only local, but it suffices for our goals, and the only tool its proof needs 
is Korn's inequality. 

Lemma 4. Let the space be at least three dimensional, D. and D as in the previous lemma. Define 
E = D\i2 iCL j wT) , where L 2 (D., TvT) c L 2 (C1, TvT) is the space oftraceless tensors. There exists a 
linear map S : L 2 (0, T) L 2 (Q, TvT) such that dom S = L 2 (£2, T), ran S c dom E n (ker EY, 
and E o S = id L 2 (a T) . This map is unique and continuous. Furthermore, ran S + c H\(Q., T). 

Proof. If S exists, it is unique. The Killing operator K was introduced at the beginning of the 
proof of the previous lemma. Here we are going to use the conformal Killing operator: K : 
L 2 (C1, T) w L 2 (Q, TVT), dom K = CS°(Q, T), K[{\ = -1/2 where l^g := £^g - 2/n V,f g. 
Similarly to the former lemma, K = E. 

The conformal version of Korn's inequality 15 states that if the space is at least three dimen- 
sional, there is a constant C such that 

liv#|| 2 < coi^ii 2 + 1|%|| 2 ) 

for any £ e C°°(0, T) (and hence for any £ e H l (D., T)). As in the proof of the previous lemma, 
Friedrichs' inequality and this one imply that for sufficiently small CI there is a y such that 

WiWlnaj) : = ll^ll 2 + ll^l 2 < rll%ll 2 , (53) 

holds for any £ e C^(C1, T) (hence for any | e H\(SX T)). Thus ker K = {0}. (In fact, ker K = 
{0}.) This is true for any metric, including the ones that admit conformal Killing vector fields 
because of the additional condition that the vector field vanishes on the boundary dCl. Hence K 
is invertible on ran K. Because of (53), this inverse is an L 2 -continuous linear map, so it can be 
extended to a continuous linear map on ran K . Therefore we have a continuous operator R on the 
entire L 2 (C1, TwT) which is defined on the closed subspace ran K by this extension, and on the 
orthogonal complement ker E = ran K 1 ^ by zero: 

R : L 2 (C1, TVT) -» J(Q, T) c L 2 (Q, T), RoK = id c » ( Q,7> ker R d ker E, R is continuous. 

We indicated that ran R c H\{C1). This follows from (53), which also implies that R is continuous 
even as map into H\{C1), but we need only its L 2 -continuity. 

Take the adjoint of R o K. All the adjoints are taken here with respect to the L 2 scalar product. 
If A and B are two densely defined operators so that AB is also densely defined, then usually only 
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(AB) + d B + A + holds, but if A is continuous and everywhere defined, then (AB) + = B + A + . This 
is our case, so (R o K) + = K + o R + . Being the adjoint of a continuous operator, R + is also 
continuous, and domR + = L 2 (D.,T). On the other hand, (R o K) + = (idc~(o,r)) + = id L 2 (ar) . 
Since K + = E, we finally have E o R + = id £ 2 (£1>r) . This also means that rani? + c dom£. As 
usual, (rani?" 1 ")- 1 - = ker R ++ = keri? D kerE. (R ++ = R because R is continuous, everywhere 
defined.) Thus ranii + c rani2 + c (ker E)- 1 . The map S := R + is exactly what we were looking 
for. □ 

Corollary. Let the space be at least three dimensional, 12 and E as in the above lemma. Let 
L : L 2 (C1, TvT) — > Rbe a continuous linear map, and ker E c ker L. Then there is an f e L 2 (I2) 
such that L[p] = (f,p) for any solution p to the constraint equation E[p] + 1/nVp = 0, where 
p := p\ € H l {H) and p := p - I In p g e dom E. 

Proof. Continuity of L implies that L = (t, •) with some t e L 2 (C1,T\/T). It is enough to prove 
the statement for traceless t, so we assume that t e L 2 {Cl,T V T). By the assumption on L, t is 
orthogonal to ker E. Any solution p to the constraint equation can be written as p = S[l/nVp] + q, 
where p = p\, the map S is the solution assignment found in the lemma, and q e ker E. Note that 
S + is defined everywhere on L 2 (0, TvT). We write L[p] = (t, S[l/n Vp] + q) = (t, S[l/n Vp]) = 
(S + [t], l/n Vp) for any p that satisfies the constraint equation. Let v := l/n S + [t] e Hl(Cl, T). 
Hl(Cl, T) is the completion of C~(0, T) in H l (Cl, T), so there is a sequence v k e C~(Q, T) that 
converges to v in H (C1,T). By the definition of the weak derivative, V can be recast on v*. 
Since V,v^ strongly converges in L 2 (D.) to the weak divergence of v, we have L[p] = (f, p) with 
/ = -V ; V 6 L 2 (H). □ 



Appendix D: Embeddability of spaces and the Hamiltonian constraint 

If the class of transverse traceless smooth tensor fields on a bounded domain fl c I is broad 
enough, then we can argue that a function / : O —> R vanishes if the orthogonality condition 
fafqVqij = holds for all fields q in this class. If the metric is flat on CI, and I is at least three 
dimensional, then the argument is short. For simplicity, let 12 c R n (n > 3). Let k e R n be an 
arbitrary vector, and choose other two unit vectors m and n so that the three vectors are mutually 
orthogonal to each other. Then q(x) := e kx (m <S> m - n<8> n) gives q^qtj = 2 e 2hx , where • is the 
standard scalar product in R". The space generated by such functions is dense in L 2 (C1), so / = 0. 

That was quick. We have to slow down if I is curved since there are no constant vector fields 
analogous to k, m, and n. It is plausible, but we cannot be sure that the conclusion of the previous 
paragraph was not merely due to the flatness of the space. Nonvanishing curvature is an obstacle 
to the construction of constant vector fields by parallel transport. Fortunately, the generalization of 
the above argument to curved spaces is not long. This statement seems to contradict the length of 
this appendix. The reason for this extended note is twofold: all the relevant theorems are included 
so that the reader is not referred to the literature for the precise statements, and the question is put 
into the context of embeddability of spaces so that what would be only a technical detail in our 
argument reveals a (disputably 18 ) important feature of Einstein's gravity. 

The freedom in the choice of a transverse momentum determines what n > 2 dimensional 
spaces can be thought of as spacial slices of an n + 1 dimensional vacuum spacetime. To see 
this, first note that if the spacetime metric solves Einstein's equations, then in particular all the 
gravitational constraints are satisfied. In terms of the extrinsic curvature K, its trace K, and the 
Ricci scalar R of the spatial metric g, the momentum and the Hamiltonian constraints are the 
equations V ( (K^ - Kg ij ) = and K ij Ky - K 2 - R = 0, respectively, or in terms of the momentum 
p = K - Kg and its trace p, Vjp ij = and p'jpy - ^ p 2 - R = 0. Actually, the solubility of 
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the constraints is also sufficient for the existence of such an embedding since if (I, g) is a smooth 
three dimensional Riemannian manifold with a smooth symmetric tensor field K on it so that all 
the constraints are satisfied, then there is a globally hyperbolic Ricci-flat spacetime with a Cauchy 
surface whose induced metric and extrinsic curvature are g and K. 19 So the real question is if we 
have enough freedom to prescribe p ij py - ^ p 2 so that the Hamiltonian constraint is satisfied by 
a transverse momentum for a given spatial metric. The Campbell-Magaard theorem asserts that 
this can be done in some neighborhood of any point of L if the latter is an analytical Riemannian 
manifold, so such n dimensional manifolds can always be locally embedded into a Ricci flat n + 1 
dimensional Riemannian or Lorentzian manifold. Of course this theorem does not provide us with 
a unique embedding, and the way in which the space is embedded into the ambient spacetime can 
be further specified. One of the most natural questions is that if the space can be embedded so that it 
is a maximal slice of the surrounding spacetime, meaning that p = on the embedded space. This 
question has been raised and answered in the affirmative in [20], however, L itself was required to 
be of Lorentzian signature. Note that this is a necessary condition, since the Hamiltonian constraint 
reads as p ij pij - R = 0, in which the first term is a positive definite expression of the traceless p 
if (I, g) is Riemannian, so the constraint cannot be satisfied for a metric whose Ricci scalar is not 
nonnegative. (The tilde over p is just our usual mnemonic for tracelessness.) 

In the arguments showing that the constraints are satisfiable for a given spatial metric, some of 
the components of the extrinsic curvature are eliminated through the algebraic equations given by 
the Hamiltonian constraint (and p = 0) so that the momentum constraint becomes an equation to 
which the Cauchy-Kowalevski theorem applies. The reason why analiticity is required and only 
local results are obtained is that the proofs rely on this theorem. We give a brief summary of 
the argument in [20] for a Riemannian space L. For later purposes we will not set the trace of 
the momentum to zero, but instead it is assumed to be prescribed by an analytic function t. We 
also want to apply the analysis to a slightly more general kinetic term Gijup ij p kl , where Gp/ = 
« (gikgji + gugjk) ~ A gijgu with constants k > and A. So the constraint quadratic in p has the 
form 2Kp ij Pij - Ap 2 = s, where s is a prescribed function. In a neighborhood of a point x e L 
we introduce Gaussian normal coordinates corresponding to a the hypersurface on which say the 
first coordinate vanishes. Note that if the metric is analytic in a local coordinate system at x, then 
there are Gaussian normal coordinates in which the metric is still analytic. In these coordinates 
the metric components gy are zero for j = 2,...,n, and gn = 1. We set the components p ij to 
zero except for p 22 , p 33 , and p lj if j = l,...,n. For this arrangement n > 3 is necessary. Together 
with p\ = p n + grip 22 + g33P 33 = (we have n + 2 constraint equations for the n + 2 nonzero 
components of pv, so we have no more freedom to specify more components of the momentum. 
After the elimination of p 33 through p\ = t, the quadratic constraint gives the following second 
order equation for p 22 : 

(54) 

where g [2 , 3] := g 22 g33 - g\y Since j =2> ...,„ is positive definite, g 22 > 0, g [2 , 3] > 0, and thus 
g 33 > 0. If the discriminant of this equation is positive, then using one branch of solutions for p 22 , 
we get the desired form for the momentum constraint: 

-TTT = 2Z2Z4(z)^-+b\z), z = (x\...,x",p n ,...,p ln ), j=l,...,n, (55) 
ax v i=ik=\ ox' 

where d[{z) and b j (z) are analytic functions on a domain that will be described soon. 

A function / is analytical on an open set £1 c R m if and only if for any compact set K c CI 
there are positive constants M and r such that / e C M , r ( x ) f° r an y x e K. f e C M , r ( x ) means that / 
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is smooth in a neighborhood of x, and 

|<9j/(x)| <Mj! r" ljl for all j e N m , 

where j is a multiindex, j! := j\ \ . . . j m \, and |j| := ji + ■ ■ ■ + j m . Now we are ready to quote the 
relevant version of the Cauchy-Kowalevski theorem. The theorem is often stated without making 
it explicit how the size of the neighborhood in which the solution must exist depends on the given 
data. For later purposes, we included such a specification of the radius of convergence, which of 
course can be seen from the classic proofs of the theorem. 21 

Theorem. (Cauchy-Kowalevski) Let a\ and b j be real analytic functions at the origin ofW +N ~ l . 
Then the system of differential equations 

Q u i n N Qyk 

-^-j = Ha{{z)— i +b J {z), z = (x 2 ,...,x n ,u u ...,u N ), j=l,...,N 

OX i=2 k= 1 OX 

with initial conditions u-> = at x l = (j = 1, . . . , N) has a system of real analytic solutions 
u\x l ,x 2 , . . . x n ) in the ball B p (0), where p depends on n, N, and on the class C M ,r(0) to which all 
the coefficients a' } k and V belong. This is the only solution that is real analytic at the origin. 

Equation (55) has a form slightly different from the expression in this theorem: the coefficients 
are allowed to depend on x l . However, this is not more general than the system in the theorem 
since one can always extend the set of unknown variables (u 1 ,..., u N ) by a new function v with a 
further equation dv/dx l = 1 and initial value v = at x l =0. Then v plays the role of x 1 . 

Let x = (0, x 2 , . . . , x n ). If s(x) is positive, and t 2 (x) is sufficiently small, then the discriminant of 
equation (54) is positive at x if p 1J (x) = 0. So the constraint equations can be reduced to the form 
(55). In order to calculate to what class C M , r (z) the coefficient functions d[ and b j in (55) belong, 
we have to evaluate the derivatives of the solution of (54) with respect to (x 1 , . . . , x'\ p 11 ,.. .,p ln ) 
at (0, x 2 , . . . , x'\ 0, . . . , 0). The derivatives with respect to (p n , . . . , p ln ) result in negative powers 
of the discriminant of (54). Combining this observation with the Cauchy-Kowalevski theorem, we 
arrive at the following characterization of local solutions to the constraints: 

Lemma 5. Let (L, g) be an at least three dimensional Riemannian manifold, (U,(f) a coordinate 
chart in which g is analytic; M, r, a, and b positive constants. S and T are two families of functions 
on U such that S, T c C\ yeU C M , r (y) in the given local coordinates; and a < s, \t\ < b for all s e S , 
t 6 T. Ifbis sufficiently small, then there is a neighborhood N of any point of U such that for all 
s e S and t e T there is a smooth symmetric tensor field p which satisfies Vtp'-i = 0, p\ = t, and 
p'ipij - crp 2 = s everywhere on N, where cr is a constant. 

This lemma guarantees a class of local solutions to our constraints on the momentum. The fol- 
lowing lemma is about the consequences of orthogonality conditions similar to the one mentioned 
at the beginning of this appendix. 

Lemma 6. Let (Z, g) be an n > 3 dimensional Riemannian manifold, x e L, (Q, (p) a coordinate 
chart around x in which g and the positive function s : O —* R. + in part b) are analytic. For any 
q > there is a bounded domain iVcfl containing x with the following property. Let f 6 L 2 (N). 
Define F : W — > R by F\ v ^ := yfgf o <p~ l and Fl^n^^ = 0. Let F be the Fourier transform ofF. 

a ) If $NfP^Pij - for any smooth transverse traceless symmetric tensor field p on N, then 
supp F n 5 e (0) = 0, where B e (0) is the ball of radius g centered at the origin. 

b) If $n fP = any momentum p that is transverse on N and satisfies p ij pij - crp 2 = s, where 
cr is a constant, then supp F Pi 5^(0) = 0. 
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Proof. Let g : R" —> R be the inverse Fourier transform of the characteristic function of B e (0), 
which is the ball of radius q centered at the origin. Such a g is analytic, so for any compact set K 
there are constants M and r such that g e C M , r iy) for any y e K. Since g is bounded, there is a 
constant c such that c + g is positive. Let g a be the translation of g by a e R n , so g a (y) = g(y - a). 

For any 6 > there is a neighborhood [/ of x such that £/, T := { }, and the family that is 
defined by S := { d } U { g a + c | |a| < 5 } in the given local coordinates satisfy the conditions of 
lemma 5. Here d can be any positive constant. Thus there is a bounded domain N <z U which has 
the property that for any s e S there is a smooth transverse traceless symmetric tensor field p that 
satisfies p'-ipy = s on N. 

Now let / be as in part a). f K „F s = for all s e S , and since S contains a constant, this implies 
that F is orthogonal to all the translations of g by \a\ < 6. Therefore 

/ die g(k)*F(k) e ika = J &"k F{k) e ika = for all \a\ < 6, 

R" \k\<e 

where • denotes the standard scalar product in R n . The function a \-> fa d ,1 kF(k) e lha is analytic 
on R" since it is the Fourier transform of a compactly supported square integrable function. Thus it 
is identically zero because the above equality shows that it vanishes on an open subset of R". That 
means that Jj"^ d"kF(k) e lk ' a = for any a e R", hence F(k) = for any \k\ < q. 

Part b) goes along the same lines. Here T := { eg a \ \a\ < 8 } (6 > 0) and S = { s } satisfy the 
conditions of lemma 5 together with some neighborhood U of x, provided that e > is sufficiently 
small. □ 



Appendix E: Proof of theorem 1 

We will use the tensor field p instead of the momentum n = sjg p. The momentum dependence is 
indicated explicitly, but the metric dependence is suppressed in our notations. The starting point is 

C[a,/3,p]-Jf?[p] + C[a,/3,p];^[p] = f / (/ |o-./>V 7 . (56) 

i 

Let p p a := pp + crq + g on a bounded domain D. D supp a U supp /3, where p and <x are constants, 
p and the traceless q are smooth symmetric tensor fields. Note that the inverse metric satisfies the 
momentum constraint: V t g ij = 0. In order for p to be a momentum, we can demand for example 
that on the entire space I it is compactly supported. Recall that we use the word momentum for 
only kinematically allowed momentum (see Appendix A), so any momentum is assumed to be 
smooth. Since the right hand side of (56) is linear in p and cr, if C'[a,/3, p pcr \-Jtf[p pcr \, which is 

p / C[a,{3, p p J Vjqj +o-f C[a,{3, p p J V jP j, (57) 
i i 

vanishes, so does the term in C[a,/3, p pa .]-Jf[p pcr ] proportional to per. The latter term is 

(nA -2k) fp C[a,0, q]-2 K f C[a,/3, g] p%, (58) 

i i 

where p = p\ and p = p - l/n p g. 

We use the operators D and E introduced in lemma 3 and 4. They map square integrable 
tensor fields on a bounded domain D. into the space of square integrable vector fields. Any limit 
or convergence will be meant in these spaces. D and E axe the closures of D and E. Recall 
that p e ker D means that there is a sequence of p k which are smooth on CI, lim^oo p k = p, and 
lim,j^oo D[p k ] = 0. Since the kernel of a closed operator is closed, ker D c ker D , but we did not 
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show that ker D = ker £),so we have to be a little careful later. Ape ker D is not necessarily 
smooth, so C or C might not be defined on it. Nevertheless, (57) is zero not only for p and q such 
that p\d 6 ker D and q\ n e ker E, but also if p\& e kerZ) and q\ n e kerE, or if p\ n e ker D and 
q\a 6 ker.E, where the evaluation of (57) and (58) on a tensor field p e ker E c ker D is meant 
by taking the limit of the integrals along an approximating sequence of momenta p k . On the other 
hand, it is unpredictable what would happen if we tried p\ n e ker D and q\ n e ker E. In this case 
(57) might not converge along momenta tending to p and q even if there is nothing pathological 
about C For instance, if they contain higher derivatives of p with smooth coefficient functions, 
there is already enough room for an unpleasant behavior. 

First let p = and p = q. If q is a smooth tensor field which is transverse and traceless on 
CI, then the first term in (58) is absent, and we have J L C[a,f3, g] q ij qij = 0. Now we show that 
this implies that C[a,/3, g] = 0. Since supp a and supp/? are compact, by using an appropriate 
partition of unity, a and can always be written as a (finite) sum of compactly supported smooth 
functions a r and /3 S such that supp a,. U supp/? s D supp C[a r ,fi s , g] is covered by one or two 
disjoint coordinate domains for any r and s, depending on how far the supports of a r and /3 S are 
from each other. So we can assume that a and /? already have this property. Let (CI, (p) be one 
of these covering charts, and K := supp a U supp/3 n CI. Assume that g is analytic in the given 
local coordinates. By part a) of lemma 6, for any q > there is an open covering of K, and thus 
by the compactness of K a finite subcovering which consists of neighborhoods N so that q and N 
have the properties as in part a ) of the lemma. By the aid of a partition of unity subordinate to the 
latter subcovering, we decompose a and /3 into the sum of a r and fi s . f L C[a r ,/3 S , g] q ij qij = for 
any smooth q that is traceless and transverse on Q. (and zero on supp a U supp/3 \ O). Part a) of 
lemma 6 with / = C[a r ,/3 S , g], the linearity of the Fourier transformation, and the arbitrariness of 
q imply that C[a,j3, g] = 0. 

So the second term in (58) is zero, and the condition on C is that limjt_»oo fz P C[a,fi, qu\ = if p 
is a momentum transverse on K, and q k is a sequence of traceless momenta such that q k \a converges 
and lim^oo E[q k \ n ] = 0. As in the former paragraph, the chart (Q, cp) is one of the (at most two) 
coordinate charts that cover supp a U supp/?, and g is assumed to be analytic in (CI, ip). (On the 
other chart, if there is any, p and q k are set to zero.) CI is chosen so that it satisfies the conditions 
of lemma 3 and 4. Let p be analytic. By part b) of lemma 6 (actually, there is no quadratic 
constraint on p here) and a decomposition of a and fi similar to the one applied in the former 
paragraph, we conclude that the above limit lim^oo f L p C[a,{3, q k ] vanishes for all analytic p. The 
contribution of C'[a,/3,p] • Jtf[p] to the right hand side of (56) is J L C'[a,/3,0] Vjpj, which also 
goes to zero along p = q k . As we have just proved, so does the contribution of C[a,/3,p] ■ Jt?[p\, 
which is f L C[a,p,p]G[g], if G[g] is analytic in (CI, <p). L aJj [p] := f a p ij = f K ty[a,fi] p ij is a 

continuous linear functional on the square integrable fields on CI. The conclusion of this paragraph 
is that ker E c ker L a ^ if G[g] is analytic in (CI, cp). 

Now we show that if p e ker D and p e H l (Cl), then p e ker L a g. The coordinate chart (CI, cp) 
has the properties as before, in particular, it satisfies the conditions of lemma 3 and 4, and as before 
we define K := supp a U supp /3 n CI. By the corollary of lemma 4, there is an f a $ e L 2 (C1) such 
that L a jj[p] = (f a j3,p) if P £ kerZ) and p e ^(CX). According to lemma 4, there is a p e ker D 
for any p 6 ^(Cl). Since ^(Cl) is dense in L 2 (C1), f a $ is unique. Note that (a,fi) i-» f a ^ is 
bilinear. What we showed in the former paragraph implies that f a p = almost everywhere on 
Cl\ K since the domain CI in the previous paragraph could be replaced by any smaller domain 
containing K. Assume that G[g] is analytic in (CI, ip), and if necessary, let us add a constant c to 
it so that it is negative on CI. Here we assume that k > 0. (k < is the same, but then G[g] + c 
should be positive.) If the momentum p is transverse on CI, and iKp^pij - Ap 2 = —j=G[g~\, then all 
the constraints J4? M [p] = are satisfied on CI, so (f a ,p,p) = 0. The next step is familiar. By part b) 
of lemma 6 and an appropriate decomposition of a and /3, we conclude that f a ^ = 0. Since there is 
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some freedom in the constant added to G[g] in order to make it negative, and f aJ} depends linearly 
on this constant, we reached the conclusion f a g = with the original analytic G[g]. 

In order to clear the way for lemma 3, we have to argue that ker D c ker L a p. This statement 
is stronger than what we proved in the former paragraph since it is imaginable that an element of 
ker D can be approximated by only such sequences p k for which the derivative Wp k of the trace and 
the divergence V,-pjjJ of the traceless part are divergent, but the sum ^ ip\ + l/nV j p k tends to zero. 
For any p e kerZ) there is a convergent sequence of momenta p k for which lim^o, D[p k ] = 0. 
Write the traceless part p k of p k as the sum of a term in ker E and a term orthogonal to ker E. Since 
ker E c ker L a ^, if p k is replaced by the latter term, L aJ3 [p k ] remains the same as with the original 
p k . So let us do this replacement. Let r k be a traceless tensor that solves E[r k ] + \/nVp k = 0, 
where p k is the trace of p k . Lemma 4 guarantees that f k exists. Define r k := f k + 1/n p k g. As 
shown in the previous paragraph, L a/3 [r k ] = 0. By the continuity of the solution assignment in 
lemma 4, lim,,^ r k = lim^oo p k since p k solves the same equation as r k apart from a term which 
converges to zero as n — » oo. Since L a g is a continuous functional, L a p[p\ = lim,,^ L at p[r k ] = 0. 
Apply lemma 3 to complete the proof, at least for a metric analytic in (O, (p). 

By properties (iii),(iv), and lemma 1, we can extend the results to any kinematically allowed 
metric. 



Appendix F 

Lemma 7. If Gfj is an ultralocal concomitant of the metric with symmetry properties Gfj = Gfj = 
Gfj, then Gfj = K(6 k i6j + d'idj) - /igijg kl , where k and A are constant. 

Proof. Take any x h-> y(x) diffeomorphism such that y(x ) = x and (y*gy)(xo) = gij(x ), that is, 
Dj := djy'(xo) e 0(n). The transformation rule of G, whose components are Gfj = Gij mn g km g ln , is 

<y.G)« r f (xo) = (D-^iD-yDttfGfjixo). 

Since G is concomitant of the metric (and ultralocal in it), (y*G(g))(xo) = G(y t g(x )). But by our 
selection of y, this implies that G(x ) is invariant, so from the transformation rule of G we get 

rf i D f .G* e f (xo) = G? j (xo)DtD[, 

which means that G(x ) is an intertwiner between representations of the orthogonal group on sym- 
metric matrices. Since the symmetrized product of two fundamental representations of 0(n) de- 
composes into two inequivalent irreducible representations, Schur's lemma tells us that the space 
of the intertwiners is only two dimensional, and the formula in the statement already spans a two 
dimensional space at every point. The coefficients k and A are constant since they are scalars, con- 
comitant of the metric, and the only way that £ l djK(xo) = whenever V ( ,-£y)(xo) = is that d t K = 0, 
and similarly, 5,vl = 0. □ 

Lemma 8. Let V be the manifold of the kinematically allowed fields on L, and G : V —* C°°(L) 
such that 

^N-G[0]=/((5 j N) Lj | <7 ), (59) 
o<p 

where f is a linear function ofNe C™(L) and its derivatives up to order J < oo with smooth 
coefficients. Then G[<p\(x) is the function of(/)(x) and the derivatives of<patx up to order J. 

Proof. Let ((f> A )Ae(o,\) 6 F be a variation. If we multiply the left hand side of (59) by d^, and 
integrate it over I, we get ^(N- G[(f> A ]). The right hand side of (59) shows that d A (N- G[<f> A ]) = 
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if supp d A <f> H supp N = since /((5jN(x))j| < y) = if x £ supp N. It can be shown that by the 
assumptions made on V - see property (i) in Appendix A - this implies that N- G[0] depends on 
restricted to any neighborhood that contains supp N. Differentiate N- G[<f> A ] with respect to A, use 
(59) for the functional derivative, integrate by parts to get undifferentiated N only, then integrate 
with respect to A as well, exchange the order of the latter integral with the integration over Z, and 
finally note that N is an arbitrary C~(Z) function. The result is 

- G[<f> Q ](x) = [ dA V 0[(f> A ](x) dfiMx) (60) 

for any point x e L. The multiindex j labels coordinates on Z. 

Let if/ e F. By property (i), there is a variation {(f> A )Ae{Q,\) e T such that (f> A = Aif/ + (1 - A)if/j on a 
neighborhood of x, and ^=^on Z. Here if/j is the Taylor series of if/ about x, truncated at order 
J. The right hand side of (60) vanishes for <f> A . So G[iffj](x) = G[<f> Q ](x) = G[(f> A=l ](x) = G[iff](x). 
G[(f> A ](x) depends only on <f> A restricted to an arbitrary neighborhood of x, and in a neighborhood 
of x the function if/j is determined by if/ and its derivatives at x up to order J. □ 
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